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Abstract 

These lecture notes provide an (almost) self-contained account on conformal 
invariance of the planar critical Ising and FK-Ising models. They present the 
theory of discrete holomorphic functions and its applications to planar statistical 
physics (more precisely to the convergence of fermionic observables) . Convergence 
to SLE is discussed briefly. Many open questions are included. 
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1 Introduction 

The celebrated Lenz-Ising model is one of the simplest models of statistical physics 
exhibiting an order-disorder transition. It was introduced by Lenz in |Len20| as an 
attempt to explain Curie's temperature for ferromagnets. In the model, iron is modeled 
as a collection of atoms with fixed positions on a crystalline lattice. Each atom has a 
magnetic spin, pointing in one of two possible directions. We will set the spin to be equal 
to 1 or -1. Each configuration of spins has an intrinsic energy, which takes into account 
the fact that neighboring sites prefer to be aligned (meaning that they have the same 
spin), exactly like magnets tend to attract or repel each other. Fix a box A c Z'^ of size 
n. Let a € {-1, 1}^ be a configuration of spins 1 or -1. The energy of the configuration 
(T is given by the Hamiltonian 

EA{cr) := - ^ a^ay 

x~y 

where x ~ y means that x and y are neighbors in A. The energy is, up to an additive 
constant, twice the number of disagreeing neighbors. Following a fundamental principle 
of physics, the spin-configuration is sampled proportionally to its Boltzmann weight: at 
an inverse-temperature f3, the probability /i/3,A of a configuration a satisfies 

e-/3SA(<T) 

where 

5-6{-14}A 

is the so-called partition function defined in such a way that the sum of the weights over 
all possible configurations equals 1. Above a certain critical inverse-temperature Pc, 
the model has a spontaneous magnetization while below /3c does not (this phenomenon 
will be described in more detail in the next section). When /3c lies strictly between 
and oo, the Ising model is said to undergo a phase transition between an ordered and 
a disordered phase. The fundamental question is to study the phase transition between 
the two regimes. 

Lenz's student Ising proved the absence of phase transition in dimension one (mean- 
ing /3c = oo) in his PhD thesis |Isi25| . wrongly conjecturing the same picture in higher 
dimensions. This belief was widely shared, and motivated Heisenberg to introduce his 
famous model |Hei28| . However, some years later Peierls |Pei36| used estimates on the 
length of interfaces between spin clusters to disprove the conjecture, showing a phase 
transition in the two-dimensional case. Later, Kramers and Wannier |KW41a[ IKW41b] 
derived nonrigorously the value of the critical temperature. 
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Figure 1: Ising configurations at (3 < f3c, at /3 = f3c, and (3> Pc respectively. 



In 1944, Onsager |Ons44j computed the partition function of the model, followed by 
further computations with Kaufman, see [KO50J for instanc(Q In the physical approach 
to statistical models, the computation of the partition function is the first step towards 
a deep understanding of the model, enabling for instance the computation of the free 
energy. The formula provided by Onsager led to an explosion in the number of results 
on the 2D Ising model (papers published on the Ising model can now be counted in the 
thousands). Among the most noteworthy results, Yang derived rigorously the sponta- 
neous magnetization |Yan52j (the result was derived nonrigorously by Onsager himself) . 
McCoy and Wu |MW73| computed many important quantities of the Ising model, includ- 
ing several critical exponents, culminating with the derivation of two-point correlations 
between sites (0, 0) and {n,n) in the whole plane. See the more recent book of Palmer 
for an exposition of these and other results |Pal07| . 

The computation of the partition function was accomplished later by several other 
methods and the model became the most prominent example of an exactly solvable 
model. The most classical techniques include the transfer-matrices technique devel- 
oped by Lieb and Baxter [ Lie67| lBax89| . the Pfaffian method, initiated by Fisher and 
Kasteleyn, using a connection with dimer models [Fis66, Kas61| . and the combinatorial 
approach to the Ising model, initiated by Kac and Ward |KW52| and then developed by 
Sherman |She60| and Vdovichenko |Vdo65| ; see also the more recent 

Despite the number of results that can be obtained using the partition function, 
the impossibility of computing it explicitly enough in finite volume made the geomet- 
ric study of the model very hard to perform while using the classical methods. The 
lack of understanding of the geometric nature of the model remained mathematically 
unsatisfying for years. 

The arrival of the renormalization group formalism (see [ Fis98| for a historical exposi- 
tion) led to a better physical and geometrical understanding, albeit mostly non-rigorous. 
It suggests that the block-spin renormalization transformation (coarse-graining, e.g. re- 
placing a block of neighboring sites by one site having a spin equal to the dominant 
spin in the block) corresponds to appropriately changing the scale and the temperature 

^This result represented a shock for the community: it was the first mathematical evidence that the 
mean-field behavior was inaccurate in low dimensions. 
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of the model. The Kramers- Wannier critical point then arises as the fixed point of the 
renormaHzation transformations. In particular, under simple rescaling the Ising model 
at the critical temperature should converge to a scaling limit, a continuous version of 
the originally discrete Ising model, corresponding to a quantum field theory. This leads 
to the idea of universality: the Ising models on different regular lattices or even more 
general planar graphs belong to the same renormaHzation space, with a unique critical 
point, and so at criticality the scaling limit and the scaling dimensions of the Ising model 
should always be the same (it should be independent of the lattice whereas the critical 
temperature depends on it). 

Being unique, the scaling limit at the critical point must satisfy translation, rotation 
and scale invariance, which allows one to deduce some information about correlations 
|PP66llKaH66] . In seminal papers |BPZ84bl lBPZ84a| . Belavin, Polyakov and Zamolod- 
chikov suggested a much stronger invariance of the model. Since the scaling-limit quan- 
tum field theory is a local field, it should be invariant by any map which is locally a 
composition of translation, rotation and homothety. Thus it becomes natural to pos- 
tulate full conformal invariance (under all conformal transformation^ of subregions) . 
This prediction generated an explosion of activity in conformal field theory, allowing 
nonrigorous explanations of many phenomena; see |ISZ88| for a collection of the original 
papers of the subject. 

To summarize, Conformal Field Theory asserts that the Ising model admits a scaling 
limit at criticality, and that this scaling limit is a conformally invariant object. From 
a mathematical perspective, this notion of conformal invariance of a model is ill-posed, 
since the meaning of scaling limit is not even clear. The following solution to this 
problem can be implemented: the scaling limit of the model could simply retain the 
information given by interfaces only. There is no reason why all the information of a 
model should be encoded into information on interfaces, yet one can hope that most of 
the relevant quantities can be recovered from it. The advantage of this approach is that 
there exists a mathematical setting for families of continuous curves. 

In the Ising model, there is a canonical way to isolate macroscopic interfaces. Con- 
sider a simply-connected domain O with two points a and b on the boundary and ap- 
proximate it by a discrete graph Us c JZ^. The boundary of $7^ determines two arcs 
dab and dba and we can fix the spins to be +1 on the arc dab and -1 on the arc dba 
(this is called Dobrushin boundary conditions). In this case, there exists an interfac^ 
separating +1 and -1 going from a to 6 and the prediction of Conformal Field Theory 
then translates into the following predictions for models: interfaces in fls converge when 
6 goes to to a random continuous non-selfcrossing curve 7(o,a,b) between a and b in D, 
which is conformally invariant in the following way: 

For any (i7,a,6) and any conformal map ip : D, ^ C, the random curve ijj o 7(57 ^,6) 
has the same law as 7(^(n),i/'(a),i/'(6))- 

In 1999, Schramm proposed a natural candidate for the possible conformally invariant 
families of continuous non-selfcrossing curves. He noticed that interfaces of models 
further satisfy the domain Markov property, which, together with the assumption of 
conformal invariance, determine the possible families of curves. In [SchOO| . he introduced 
the Schramm- Loewner Evolution (SLE for short): for k > 0, the SLE(k;) is the random 
Loewner Evolution with driving process ^/nBt, where (-Bt) is a standard Brownian 
motion (see Beffara's course in this volume). In our case, it implies that the random 
continuous curve 7(n,a,6) described previously should be an SLE. 

Proving convergence of interfaces to an SLE is fundamental. Indeed, SLE processes 
are now well-understood and their path properties can be related to fractal properties 

^j. e. one-to-one holomorphic maps. 

^In fact the interface is not unique. In order to solve this issue, consider the closest interface to d^b- 
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Figure 2: An interface between + and - in the Ising model. 



of the critical phase. Critical exponents can then be deduced from these properties via 
the so-called scaling relations. These notes provide an (almost) self-contained proof of 
convergence to SLE for the two-dimensional Ising model and its random-cluster repre- 
sentation the FK-Ising model (see Section |3] for a formal definition) . 



Main result 1 (Theorem 2.10) The law of interfaces of the critical Ising model converges 
in the scaling limit to a conformally invariant limit described by the Schramm- Loewner 
Evolution of parameter k = 3. 



Main result 2 (Theorem 3.13) The law of interfaces of the critical FK-Ising model 
converges in the scaling limit to a conformally invariant limit described by the Schramm- 
Loewner Evolution of parameter k = 16/3. 

Even though we now have a mathematical framework for conformal invariance, it 
remains difficult to prove convergence of interfaces to SLEs. Observe that working with 
interfaces offers a further simplification: properties of these interfaces should also be 
conformally invariant. Therefore, one could simply look at a discrete observable of the 
model and try to prove that it converges in the scaling limit to a conformally covariant 
object. Of course, it is not clear that this observable would tell us anything about critical 
exponents, yet it already represents a significant step toward conformal invariance. 

In 1994, Langlands, Pouliot and Saint- Aubin |LPSA94] published a number of numer- 
ical values in favor of conformal invariance (in the scaling limit) of crossing probabilities 
in the percolation model. More precisely, they checked that, taking different topological 
rectangles, the probability Cs{^, A, B,C, D) of having a path of adjacent open edges 
from AB to CD converges when S goes to towards a limit which is the same for 
{il, A, B,C, D) and {il' , A' , B' ,C' , D') if they are images of each other by a conformal 
map. The paper [LPSA94) . while only numerical, attracted many mathematicians to 
the domain. The same year, Cardy |Car92| proposed an explicit formula for the limit of 
percolation crossing probabilities. In 2001, Smirnov proved Cardy's formula rigorously 
for critical site percolation on the triangular lattice |Smi01| . hence rigorously provid- 
ing a concrete example of a conformally invariant property of the model. A somewhat 
incredible consequence of this theorem is that the mechanism can be reversed: even 
though Cardy's formula seems much weaker than convergence to SLE, they are actually 
equivalent. In other words, conformal covariance of one well-chosen observable of the 
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model can be sufficient to prove conformal invariance of interfaces. 

It is also possible to find an observable with this property in the Ising case (see 
Definition 2.9 1. This observable, called the fermionic observable, is defined in terms 
of the so-called high temperature expansion of the Ising model. Specific combinatorial 
properties of the Ising model translate into local relations for the fermionic observable. 
In particular, the observable can be proved to converge when taking the scaling limit. 
This convergence result (Theorem 2.111 is the main step in the proof of conformal 
invariance. Similarly, a fermionic observable can be defined in the FK-Ising case, and 
its convergence implies the convergence of interfaces. 

Archetypical examples of conformally covariant objects are holomorphic solutions to 
boundary value problems such as Dirichlet or Riemann problems. It becomes natural 
to expect that discrete observables which are conformally covariant in the scaling limit 
are naturally preharmonic or preholomorphic functions, i.e. relevant discretizations of 
harmonic and holomorphic functions. Therefore, the proofs of conformal invariance har- 
ness discrete complex analysis in a substantial way. The use of discrete holomorphicity 
appeared first in the case of dimers |KenOO| and has been extended to several statistical 
physics models since then. Other than being interesting in themselves, preholomorphic 
functions have found several applications in geometry, analysis, combinatorics, and prob- 
ability. We refer the interested reader to the expositions by Lovasz |Lov04| , Stephenson 
|Ste05| . Mercat |Mer01| . Bobenko and Suris f BS08| . Let us finish by mentioning that the 
previous discussion sheds a new light on both approaches described above: combinato- 
rial properties of the discrete Ising model allow us to prove the convergence of discrete 
observables to conformally covariant objects. In other words, exact integrability and 
Conformal Field Theory are connected via the proof of the conformal invariance of the 
Ising model. 
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1.1 Organization of the notes 

Section|2]presents the necessary background on the spin Ising model. In the first subsec- 
tion, we recall general facts on the Ising model. In the second subsection, we introduce 
the low and high temperature expansions, as well as Kramers- Wannier duality. In the 
last subsection, we use the high-temperature expansion in spin Dobrushin domains to 
define the spin fermionic observable. Via the Kramers- Wannier duality, we explain how 
it relates to interfaces of the Ising model at criticality and we state the conformal in- 
variance result for Ising. 

Section[3]introduces the FK-Ising model. We start by defining general FK percolation 
models and we discuss planar duality. Then, we explain the Edwards-Sokal coupling, an 
important tool relating the spin Ising and FK-Ising models. Finally, we introduce the 
loop representation of the FK-Ising model in FK Dobrushin domains. It allows us to 



6 



define the FK fermionic observable and to state the confornial invariance result for the 
FK-Ising model. 

Section |4] is a brief survey of discrete complex analysis. We first deal with prehar- 
monic functions and a few of their elementary properties. These properties will be used 
in Section |6] In the second subsection, we present a brief historic of preholomorphic 
functions. The third subsection is the most important, it contains the definition and 
several properties of s-holomorphic (or spin-holomorphic) functions. This notion is cru- 
cial in the proof of conformal invariance: the fermionic observables will be proved to 
be s-holomorphic, a fact which implies their convergence in the scaling limit. We also 
include a brief discussion on complex analysis on general graphs. 

Section |5] is devoted to the convergence of the fermionic observables. First, we show 
that the FK fermionic observable is s-holomorphic and that it converges in the scaling 
limit. Second, we deal with the spin fermionic observable. We prove its s-holomorphicity 
and sketch the proof of its convergence. 

Section [6] shows how to harness the convergence of fermionic observables in order 
to prove conformal invariance of interfaces in the spin and FK-Ising models. It mostly 
relies on tightness results and certain properties of Loewner chains. 

Section [7] is intended to present several other applications of the fermionic ob- 
servables. In particular, we provide an elementary derivation of the critical inverse- 
temperature. 

Section [8] contains a discussion on generalizations of this approach to lattice models. 
It includes a subsection on the Ising model on general planar graphs. It also gathers 
conjectures regarding models more general than the Ising model. 

1.2 Notations 

1.2.1 Primal, dual and medial graphs 

We mostly consider the (rotated) square lattice L with vertex set e'''^l^l? and edges 
between nearest neighbors. An edge with end-points x and y will be denoted by [xy]. 
If there exists an edge e such that e = [xy], we write x ~ y. Finite graphs G will always 
be subgraphs of L and will be called primal graphs. The boundary of G, denoted by 
dG, will be the set of sites of G with fewer than four neighbors in G. 

The dual graph G* of a planar graph G is defined as follows: sites of G* correspond 
to faces of G (for convenience, the infinite face will not correspond to a dual site), edges 
of G* connect sites corresponding to two adjacent faces of G. The dual lattice of L is 
denoted by L*. 

The medial lattice L* is the graph with vertex set being the centers of edges of L, 
and edges connecting nearest vertices, see Fig.|6] The medial graph is the subgraph 
of composed of all the vertices of L* corresponding to edges of G. Note that is 
a rotated and rescaled (by a factor l/\/2) version of L, and that it is the usual square 
lattice. We will often use the connection between the faces of and the sites of L and 
L* . We say that a face of the medial lattice is black if it corresponds to a vertex of L, 
and white otherwise. Edges of L* are oriented counterclockwise around black faces. 

1.2.2 Approximations of domains 

We will be interested in finer and finer graphs approximating continuous domains. For 
i5 > 0, the square lattice \/2Sh of mesh-size \/25 will be denoted by . The definitions 
of dual and medial lattices extend to this context. Note that the medial lattice has 
mesh-size 6. 

For a simply connected domain 17 in the plane, we set fi^ = $7 n L^. The edges 
connecting sites of 17^ are those included in Q.. The graph fi^ should be thought of as 
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a discretization of fl (we avoid tccliriicalitics concerning the regularity of the domain). 
More generally, when no continuous domain is specified, fls stands for a finite simply 
connected (meaning that the complement is connected) subgraph of L^. 

We will be considering sequences of functions on for S going to 0. In order to 
make functions live in the same space, we implicitly perform the following operation: 
for a function / on O^, we choose for each square a diagonal and extend the function to 
f2 in a piecewise linear way on every triangle (any reasonable way would do). Since no 
confusion will be possible, we denote the extension by / as well. 



1.2.3 Distances and convergence 

Points in the plane will be denoted by their complex coordinates, Ile{z) and Im(^) will 
be the real and imaginary parts of z respectively. The norm will be the usual complex 
modulus I • |. Unless otherwise stated, distances between points (even if they belong to 
a graph) are distances in the plane. The distance between a point z and a closed set 
F is defined by 

diz,F) := inf|z-y|. (1.1) 

Convergence of random parametrized curves (say with time-parameter in [0, 1]) is in 
the sense of the weak topology inherited from the following distance on curves: 

rf(7i,72) = inf sup |7i(w) -72(</'(w))|, (1-2) 

V U€[0,1] 

where the infimum is taken over all reparametrizations (i.e. strictly increasing continu- 
ous functions (j): [0, 1] ^ [0, 1] with 0(0) = and </>(!) = 1). 



2 Two-dimensional Ising model 

2.1 Boundary conditions, infinite-volume measures and phcise 
transition 

The (spin) Ising model can be defined on any graph. However, we will restrict ourselves 
to the (rotated) square lattice. Let G be a finite subgraph of L, and b € {-1,+1}^'^. 
The Ising model with boundary conditions 6 is a random assignment of spins {-1,+1} 
(or simply -/+) to vertices of G such that cTx = b^ on dG, where cr^ denotes the spin at 
site X. The partition function of the model is denoted by 
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exp 



<7€{-l,l}G: o-=6 on dG 



L x~y 



(2.1) 



where j3 is the inverse-temperature of the model and the second summation is over all 
pairs of neighboring sites x,y in G. The probability of a configuration a is then equal 
to 



L x~y 



(2.2) 



Equivalently, one can define the Ising model without boundary conditions, also called 
free boundary conditions (it is the one defined in the introduction). The measure with 
free boundary conditions is denoted by q. 
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We will not offer a complete exposition on the Ising model and we rather focus on 
crucial properties. The following result belongs to the folklore (see |FKG71| for the 
original paper) . An event is called increasing if it is preserved by switching some spins 
from - to +. 

Theorem 2.1 (Positive association at every temperature). The Ising model on a finite 
graph G at temperature /3 > satisfies the following properties: 

• FKG inequality: For any boundary conditions h and any increasing events A, B, 

4,G{AnB) > (2.3) 

• Comparison between boundary conditions: For boundary conditions hi < &2 
(meaning that spins + in hi are also + m and an increasing event A, 

If ( |2.4| is satisfied for every increasing event, we say that /i^^^ stochastically dom- 
inates ^J■''^Q (denoted by /^J^^g ^ l^^pc^' Two boundary conditions are extremal for the 
stochastic ordering: the measure with all + (resp. all -) boundary conditions, denoted 
by ^ip Q (resp. A'^,g) largest (resp. smallest). 

Theorem 12.11 enables us to define infinite- volume measures as follows. Consider the 
nested sequence of boxes A„ = [-n,n]^. For any > and any increasing event A 
depending only on spins in Ajv, the sequence (/k^ A„(^))ra>Ar is decreasing The limit, 
denoted by /i^(A), can be defined and verified to be independent on A^. 

In this way, /i^ is defined for increasing events depending on a finite number of 
sites. It can be further extended to a probability measure on the cr-algebra spanned 
by cylindrical events (events measurable in terms of a finite number of spins). The 
resulting measure, denoted by /i^, is called the infinite- volume Ising model with + 
boundary conditions. 

Observe that one could construct (a priori) different infinite-volume measures, for 
instance with - boundary conditions (the corresponding measure is denoted by /i^). If 
infinite- volume measures are defined from a property of compatibility with finite volume 
measures, then /.t^ and /i^ are extremal among infinite- volume measures of parameter 
/3. In particular, if = /^^, there exists a unique infinite volume measure. 

The Ising model in infinite- volume exhibits a phase transition at some critical inverse- 
temperature j3c- 

Theorem 2.2. Let = |ln(l + \/2). The magnetization /^^[o'o] at the origin is strictly 
positive for j3 > [3^ and equal to when (5 < /3c- 

In other words, when /? > /3c, there is long range memory, the phase is ordered. When 
/3 < /3c, the phase is called disordered. The existence of a critical temperature separating 
the ordered from the disordered phase is a relatively easy fact pei36j (although at the 
time it was quite unexpected). Its computation is more difficult. It was identified with- 
out proof by Kramers and Wannier [ KW41a|, IKW41b ] using the duality between low and 
high temperature expansions of the Ising model (see the argument in the next section). 
The first rigorous derivation is due to Yang |Yan52| . He uses Onsager's exact formula 
for the (infinite-volume) partition function to compute the spontaneous magnetization 
of the model. This quantity provides one criterion for localizing the critical point. The 

^Indeed, for any configuration of spins in dAn being smaller than all +, the restriction of /^^ ^ to 
An is stochastically dominated by ^ . 



9 



first probabilistic computation of the critical inverse-temperature is due to Aizenman, 



Barsky and Fernandez |ABF87| . In Subsection 7.1 we present a short alternative proof 
of Theorem |2. 2 1 using the fermionic observable. 

The critical inverse-temperature has also an interpretation in terms of infinite- volume 
measures (these measures are called Gibbs measures). For /? < /3c there exists a unique 
Gibbs measure, while for /3 > /3c there exist several. The classification of Gibbs measures 
in the ordered phase is interesting: in dimension two, any infinite- volume measure is a 



convex combination of /i^ and /i^ (see [ Aiz80[ Hig81 or the recent proof |CV10| ). This 



result is no longer true in higher dimension: non-translational-invariant Gibbs measures 
can be constructed using 3D Dobrushin domains |Dob72| . 

When /? > /3c, spin-correlations ^^^[(Joax] do not go to when x goes to infinity. 
There is long range memory. At /3c, spin-correlations decay to following a power law 
|()ns44| : 

when X ^ oo. When (3 < l3c, spin-correlations decay exponentially fast in \x\. More 
precisely, we will show the following result first due to |MW73| : 

Theorem 2.3. For P < /3c, and a = e"/'*(a; + iy) e C, 

r«(a) = lim — In/igrcrocrnal] = ^ arcsinhfsa:) + y arcsinhfsy) 
where [no] is the site o/L closest to na, and s solves the equation 

^/^T^2^ + yiT^V = sinh(2/3) + (sinh(2/3))"^ . 

The quantity Tp{z) is called the correlation length in direction z. When getting 
closer to the critical point, the correlation length goes to infinity and becomes isotropic 
(it does not depend on the direction, thus giving a glimpse of rotational invariance at 
criticality) : 

Theorem 2.4 (see e.g. |Mes06| ). For z e C, the correlation length satisfies the following 
equality 

lim-J^=4N. (2.5) 

13^ Pa (/3c -/3) 

2.2 Low and high temperature expansions of the Ising model 

The low temperature expansion of the Ising model is a graphical representation on the 
dual lattice. Fix a spin configuration a for the Ising model on G with + boundary 
conditions. The collection of contours of a spin configuration a is the set of interfaces 
(edges of the dual graph) separating + and - clusters. In a collection of contours, an 
even number of dual edges automatically emanates from each dual vertex. Reciprocally, 
any family of dual edges with an even number of edges emanating from each dual vertex 
is the collection of contours of exactly one spin configuration (since we fix + boundary 
conditions) . 

The interesting feature of the low temperature expansion is that properties of the 
Ising model can be restated in terms of this graphical representation. We only give the 
example of the partition function on G but other quantities can be computed similarly. 
Let £g* be the set of possible collections of contours, and let \uj\ be the number of edges 
of a collection of contours cj, then 

^ = ^'^s''^ ^* ^ (e-2^)'"'. (2.6) 
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The high temperature expansion of the Ising model is a graphical representation on 
the primal lattice itself. It is not a geometric representation since one cannot map a spin 
configuration ct to a subset of configurations in the graphical representation, but rather 
a convenient way to represent correlations between spins using statistics of contours. It 
is based on the following identity: 

e/3<Tx<T„ ^ cosh(/3) + a^ay sinh(/3) = cosh(/3) [1 + tanh(/3)cr^crj,] (2.7) 

Proposition 2.5. Let G be a finite graph and a, b be two sites of G. At inverse- 
temperature (3 > 0, 

zi.G= 2# ^'='^"'='=^^cosh(/3)# "'scsinG ^ tanh(;3)l'^l (2.8) 

/ . . £^.go(a.6)tanh(/?)l'-l 
^^'-^'^'•^^^^ E....tanh(/3)H ' (^.9) 

where £q (resp. £Q{a,b) ) is the set of families of edges of G such that an even number 
of edges emanates from each vertex (resp. except at a and b, where an odd number of 
edges emanates). 

The notation coincides with the definition Sq* in the low temperature expansion 
for the dual lattice. 



Proof Let us start with the partition function (2.8). Let E be the set of edges of G. 
We know 



E n e^'^^^" 

cr [xy]eE 

cosh(/3)#<='is<=^''^«^ n [l + tanh(/3)a,aj,] 

cr [xy]eE 

cosh(/3)# -'iscs in G ^ ^ tanh(/3)l-l fl '^-^y 

o" uj(zE e=[xy]euj 

cosh(/3)#<='is<=^-« ^ tanh(/?)l'^l^ U a^ay 



where we used (2.7) in the second equality. Notice that Y,(7Y\e=[xy]<iui^x<^y equals 
2# vortices G •£ 1^ g^j-^^ g otherwise, hence proving (2.8). 

Fix a,b e G. By definition. 



Z,e"/3^(-) 



(2.10) 



where H{(t) = -T,i~j o'icrj. The second identity boils down to proving that the right 
hand terms of (2.9) and (2.10) are equal, i.e. 



Ea.a,e-^^('^) = 2* 



vertices G cosh(/3)# ^'^Ses in G ^ tanh(^)l' 



(2.11) 



The first lines of the computation for the partition function are the same, and we end 
up with 



-^«('^) = cosh(/3)#<='is<=^-« ^ tanh(/?)l"l^a,a6 fl <^-'^y 
= 2#™'"™"'^cosh(/3)#'='^sosinG ^ tanh(^)l"l 

a;e5G(a,t)) 
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Figure 3: The possible collections of contours for + boundary conditions in the low- 
temperature expansion do not contain edges between boundary sites of G. Therefore, 
they correspond to collections of contours in Eq*, which are exactly the collection of 
contours involved in the high-temperature expansion of the Ising model on G* with free 
boundary conditions. 



since Xlo- o'aO'b ne=[2:!;]etj <7x<^y cquals 2^ vertices G jf ^ g ^^(a, 5), and Otherwise. □ 

The set £g is the set of collections of loops on G when forgetting the way we draw 
loops (since some elements of £g, like a figure eight, can be decomposed into loops in 
several ways), while £0(0,, b) is the set of collections of loops on G together with one 
curve from a to b. 

Proposition 2.6 (Kramers- Wannier duality). Let /? > and define /?* e (0, 00) such 
that tanh(/3*) = e^'^^ , then for every graph G, 

2 * ^-^^-^^ <=* cosh(/3*) # ^"^"^^ - = {e^f " 4.,g* ■ (2-12) 

Proof When writing the contour of connected components for the Ising model with + 
boundary conditions, the only edges of L* used are those of G* . Indeed, edges between 
boundary sites cannot be present since boundary spins are +. Thus, the right and left- 



hand side terms of (2.121 both correspond to the sum on £g* of (c 2^)l'^l or equivalently 
of tanh(/3*)l"', implying the equality (see Fig.js]). □ 

We are now in a position to present the argument of Kramers and Wannier. Physicists 
expect the partition function to exhibit only one singularity, localized at the critical 
point. If f3* + j3c, there would be at least two singularities, at /3c and /?*, thanks to the 
previous relation between partition functions at these two temperatures. Thus, f3c must 
equal /?*, which implies /3c = + \/2). Of course, the assumption that there is a 

unique singularity is hard to justify. 

Exercise 2.7. Extend the low and high temperature expansions to free and + boundary 
conditions respectively. Extend the high-temperature expansion to n-point spin correla- 
tions. 

Exercise 2.8 (Peierls argument). Use the low and high temperature expansions to show 
that Pc e (0,00), and that correlations between spins decay exponentially fast when (3 is 
small enough. 
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Figure 4: An example of collection of contours in £{as, zs) on the lattice fio- 



2.3 Spin-Dobrushin domain, fermionic observable and results on 
the Ising model 

In this section we discuss the scaling limit of a single interface between + and - at 
criticality. We introduce the fundamental notions of Dobrushin domains and the so- 
called fermionic observable. 

Let (S7, a, b) be a simply connected domain with two marked points on the boundary. 
Let be the medial graph of ^Ig composed of all the vertices of bordering a black 
face associated to ^s, see Fig|4] This definition is non-standard since we include medial 
vertices not associated to edges of fls- Let and bg be two vertices of dftg close to a 
and b. We further require that bs is the southeast corner of a black face. We call the 
triplet (flg,as,bs) a spin- Dobrushin domain. 

Let Zs s ^g- Mimicking the high-temperature expansion of the Ising model on Qg, 
let £{ag,zg) be the set of collections of contours drawn on fig composed of loops and 
one interface from ag to zg, see Fig. |4] For a loop configuration uj, 7(0;) denotes the 
unique curve from ag to zg turning always left when there is an ambiguity. With these 
notations, we can define the spin-Ising fermionic observable. 

Definition 2.9. On a spin Dobrushin domain {flg,ag,bg), the spin-Ising fermionic ob- 
servable at Zg e is defined by 



where the winding W-y{ag, zg) is the (signed) total rotation in radians of the curve 7 
between ag and zg. 

The complex modulus of the denominator of the fermionic observable is connected 
to the partition function of a conditioned critical Ising model. Indeed, fix bg e d^\. 
Even though E{ag^bg) is not exactly a high-temperature expansion (since there are two 
half-edges starting from ag and bg respectively), it is in bijection with the set f (a,5). 
Therefore, ( 2.11[ ) can be used to relate the denominator of the fermionic observable 
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Figure 5: A high temperature expansion of an Ising model on the primal lattice together 
with the corresponding configuration on the dual lattice. The constraint that ag is con- 
nected to bg corresponds to the partition function of the Ising model with +/- boundary 
conditions on the domain. 



to the partition function of the Ising model on the primal graph with free boundary 
conditions conditioned on the fact that a and h have the same spin. Let us mention that 
the numerator of the observable also has an interpretation in terms of disorder operators 
of the critical Ising model. 

The weights of edges are critical (since \/2 - 1 = e"^^"^). Therefore, the Kramers- 
Wannier duality has an enlightening interpretation here. The high-temperature expan- 
sion can be thought of as the low-temperature expansion of an Ising model on the dual 
graph, where the dual graph is constructed by adding one layer of dual vertices around 
dG, see Fig. |5] Now, the existence of a curve between as and hs is equivalent to the 
existence of an interface between pluses and minuses in this new Ising model. Therefore, 
it corresponds to a model with Dobrushin boundary conditions on the dual graph. This 
fact is not surprising since the dual boundary conditions of the free boundary conditions 
conditioned on ctq = cTf, are the Dobrushin ones. 

From now on, the Ising model on a spin Dobrushin domain is the critical Ising model 
on with Dobrushin boundary conditions. The previous paragraph suggests a connec- 
tion between the fermionic observable and the interface in this model. In fact. Section [6] 
will show that the fermionic observable is crucial in the proof that the unique interface 
75 going from as to hs between the + component connected to the arc 9*^^ and the - 
component connected to (preserve the convention that the interface turns left every 
time there is a choice) is conformally invariant in the scaling limit. Figures [l] (center 
picture) and |2] show two interfaces in domains with Dobrushin boundary conditions. 

Theorem 2.10. Let {fl,a,b) be a simply connected domain with two marked points on 
the boundary. Let 75 be the interface of the critical Ising model with Dobrushin boundary 
conditions on the spin Dobrushin domain (Jl^, 05,65). Then {•^s)s>a converges weakly 
as 6 ^ to the (chordal) Schramm- Loewner Evolution with parameter k = 3. 

The proof of Theorem |2.10| follows the program below, see Section |6] 
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• Prove that the family of interfaces {■ys)s>o is tight. 

• Prove that M['^ = ^si«\^5[o,t],7a(t),&a (25) is a martingale for the discrete curve jg. 

• Prove that these martingales are converging when S goes to 0. This provides us 
with a continuous martingale {M^)t for any sub-sequential limit of the family 
{75)s>o- 

• Use the martingales {M^)t to identify the possible sub-sequential limits. Actually, 
we will prove that the (chordal) Schramm-Loewner Evolution with parameter k = 3 
is the only possible limit, thus proving the convergence. 

The third step (convergence of the observable) will be crucial for the success of this 
program. We state it as a theorem on its own. The connection with the other steps will 
be explained in detail in Section [6j 

Theorem 2.11 ( |CS09| ). Let fl be a simply connected domain and a,b two marked 
points on its boundary, assuming that the boundary is smooth in a neighborhood of b. 
We have that 



Psis,as,bs (■) 

uniformly on every compact subset of fl, where ijj is any conformal map from Q to the 
upper half-plane H, mapping a to 00 and b to 0. 

The fermionic observable is a powerful tool to prove conformal invariance, yet it 
is also interesting in itself. Being defined in terms of the high-temperature expansion 
of the Ising model, it expresses directly quantities of the model. For instance, we will 
explain in Section |6] how a more general convergence result for the observable enables 
us to compute the energy density. 

Theorem 2.12 ( pS10| ). Let ft be a simply connected domain and aeil. If es = [xy] 
denotes the edge of 5Z^ closest to a, then the following equality holds: 

^'k.nn5z^ I'^^^'^yi = ^ —6 + o{6), 

where fi^ is the Ising measure at criticality and (f>a is the unique conformal map from 
n to the disk D sending a to and such that 4>'a{a.) > 0. 

3 Two-dimensional FK-Ising model 

In this section, another graphical representation of the Ising model, called the FK-Ising 
model, is presented in detail. Its properties will be used to describe properties of the 
Ising model in the following sections. 

3.1 FK percolation 

We refer to |Gri06| for a complete study on FK percolation (invented by Fortuin and 
Kasteleyn |FK72| ). A configuration w on G is a random subgraph of G, composed of 
the same sites and a subset of its edges. The edges belonging to lu are called open, the 
others closed. Two sites x and y are said to be connected (denoted by a; <->■ y), if there 
is an open path — a path composed of open edges — connecting them. The maximal 
connected components are called clusters. 



when (5 -> 



(2.13) 
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Boundary conditions ^ are given by a partition of dG. Let o(a;) (resp. c{uj)) denote 
the number of open (resp. closed) edges of uj and the number of connected 

components of the graph obtained from w by identifying (or wiring) the vertices in ^ 
that belong to the same class of ^. 

The FK percolation (/i^ ^ on a finite graph G with parameters p e [0, 1], and q e 
(0,00) and boundary conditions ^ is defined by 



pO(c^)(l_p)c(^)gfe('^,C) 



{u) := '—^ ^ , (3.1) 



for any subgraph w of G, where ^ q is a normalizing constant called the partition 
function for the FK percolation. Here and in the following, we drop the dependence on 
£, in k{uj,0- 

The FK percolations with parameter q < I and q> 1 behave very differently. For now, 
we restrict ourselves to the second case. When q > 1, the FK percolation is positively 
correlated: an event is called increasing if it is preserved by addition of open edges. 

Theorem 3.1. For q> 1 and p € [0, 1], the FK percolation on G satisfies the following 
two properties: 

• FKG inequality: For any boundary conditions ^ and any increasing events A, B, 

<,,G(^n^) ^ <,,g(^)<,,g(^)- (3-2) 

• Comparison between boundary conditions: for any ^ refinement of tp and 
any increasing event A, 

^ 4,,,g(^)- (3-3) 



The previous result is very similar to Theorem |2.1[ As in the Ising model case, one 
can define a notion of stochastic domination. Two boundary conditions play a special 
role in the study of FK percolation: the wired boundary conditions, denoted by ^ = 1, 
are specified by the fact that all the vertices on the boundary are pairwise connected. 
The free boundary conditions, denoted by f = 0, are specified by the absence of wirings 
between boundary sites. The free and wired boundary conditions are extremal among 
all boundary conditions for stochastic ordering. 

Infinite-volume measures can be defined as limits of measures on nested boxes. In 
particular, we set (p^ ^ for the infinite- volume measure with wired boundary conditions 
and ^ for the infinite-volume measure with free boundary conditions. Like the Ising 
model, the model exhibits a phase transition in the infinite-volume limit. 

Theorem 3.2. For any q> 1, there exists Pc{q) e (0, 1) such that for any infinite volume 
measure (pp^q, 

• if p < Pcil), there is almost surely no infinite cluster under (j)p^q, 

• if p > Pcil), there is almost surely a unique infinite cluster under (j)p,q. 

Note that g = 1 is simply bond percolation. In this case, the existence of a phase 
transition is a well-known fact. The existence of a critical point in the general case q>l 
is not much harder to prove: a coupling between two measures (pp^^q^ and 4>p2,q,G can 
be constructed in such a way that (pp-^^q.c stochastically dominates 4'p2,q,G if Pi ^ P2 
(this coupling is not as straightforward as in the percolation case, see e.g. |Gri06| ). The 
determination of the critical value is a much harder task. 
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A natural notion of duality also exists for the FK percolation on the square lattice 
(and more generally on any planar graph). We present duality in the simplest case of 
wired boundary conditions. Construct a model on G* by declaring any edge of the dual 
graph to be open (resp. closed) if the corresponding edge of the primal graph is closed 
(resp. open) for the initial FK percolation model. 

Proposition 3.3. The dual model of the FK percolation with parameters {p, q) with 
wired boundary conditions is the FK percolation with parameters {p* , q) and free bound- 
ary conditions on G* , where 

P=v{p,q) == ,/^"w (3.4) 
{l-p)q + p 

Proof Note that the state of edges between two sites of dG is not relevant when 

boundary conditions are wired. Indeed, sites on the boundary are connected via bound- 
ary conditions anyway, so that the state of each boundary edge does not alter the 
connectivity properties of the subgraph, and is independent of other edges. For this 
reason, forget about edges between boundary sites and consider only inner edges (which 
correspond to edges of G*): o(w) and c{(jj) then denote the number of open and closed 
inner edges. 

Set e* for the dual edge of G* associated to the (inner) edge e. From the definition of 
the dual configuration w* of w, we have o(a;*) = a-o{ui) where a is the number of edges 
in G* and o(a;*) is the number of open dual edges. Moreover, connected components of 
w* correspond exactly to faces of a;, so that f{uj) = k{oj*), where /(w) is the number of 
faces (counting the infinite face). Using Euler's formula 

# edges + # connected components + 1 = #sites + # faces, 

which is valid for any planar graph, we obtain, with s being the number of sites in G, 

k(co) = s - 1 + f(uj) - o(u)) = s - 1 + k(co*) - a + o(oj*). 

The probability of co* is equal to the probability of co under (pQ p g, i-c. 

= ^-^^[qii-p)/pr^^*^q''^''*^ = ^i^^qM'^n 

since g(l -p)/p = P*/(l -p*), which is exactly the statement. □ 

It is then natural to define the self-dual point Psd = Psdio) solving the equation 
P*sd = Psd, which gives 

Psd = Psd{q) - — -■ 

Note that, mimicking the Kramers- Wannier argument, one can give a simple heuristic 
justification in favor oipc{q) = Psd{q)- Recently, the computation oipc{q) was performed 
for every q>l: 
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Theorem 3.4 f |BDC10] l. The critical parameter Pc{q) of the FK percolation on the 
square lattice equals Psd{q) = \/ql{^ + \/q) for every q>l. 

Exercise 3.5. Describe the dual of a FK percolation with parameters {p,q) and free 
boundary conditions. What is the dual model of the FK percolation in infinite-volume 
with wired boundary conditions? 

Exercise 3.6 (Zhang's argument for FK percolation, |Gri06p . Consider the FK per- 
colation with parameters q > 1 and p = Psdiq)- We suppose known the fact that infinite 
clusters are unique, and that the probability that there is an infinite cluster is or 1. 
Assume that there is a.s. an infinite cluster for the measure 4>^^^ q- 

1) Let £ < 1/100. Show that there exists n > such that the (f>p^^ ^-probability that 
the infinite cluster touches [-n,n]^ is larger than 1 - e. Using the FKG inequality for 
decreasing events (one can check that the FKG inequality holds for decreasing events as 
well), show that the (f>p^^ ^-probability that the infinite cluster touches {n} x [-n,n] from 
the outside of [-n, n]^ is larger than 1 - £i . 

2) Using the uniqueness of the infinite cluster and the fact that the probability that 
there exists an infinite cluster equals or 1 (can you prove these facts?), show that a.s. 
there is no infinite cluster for the FK percolation with free boundary conditions at the 
self-dual point. 

3) Is the previous result necessarily true for the FK percolation with wired boundary 
conditions at the self- dual point? What can be said about pc{q)? 

Exercise 3.7. Prove Euler's formula. 

3.2 FK-Ising model and Edwards-Sokal coupling 

The Ising model can be coupled to the FK percolation with cluster- weight q = 2 [ES88J . 
For this reason, the q = 2 FK percolation model will be called the FK-Ising model. 
We now present this coupling, called the Edwards-Sokal coupling, along with some 
consequences for the Ising model. 

Let G be a finite graph and let a; be a configuration of open and closed edges on G. 
A spin configuration a can be constructed on the graph G by assigning independently 
to each cluster of a; a + or - spin with probability 1/2 (note that all the sites of a cluster 
receive the same spin). 

Proposition 3.8. Letp e (0, 1) and G a finite graph. If the configuration lo is distributed 
according to a FK measure with parameters (p, 2) and free boundary conditions, then the 
spin configuration a is distributed according to an Ising measure with inverse-temperature 
/3 = - | ln(l - p) and free boundary conditions. 

Proof Consider a finite graph G, let p e (0, 1). Consider a measure P on pairs (w, cr), 
where a; is a FK configuration with free boundary conditions and a is the corresponding 
random spin configuration, constructed as explained above. Then, for (cjjCr), we have: 

P[{uj,ct)] = _J_p°('^)(i_p)'=('^)2'^-(").2-'^('^) = -^p°("')(l-p)"<"\ 

P,2,G P,2,G 

Now, we construct another measure P on pairs of percolation configurations and spin 
configurations as follows. Let ct be a spin configuration distributed according to an Ising 
model with inverse-temperature /3 satisfying e"^^ = 1 - p and free boundary conditions. 
We deduce ui from a by closing all edges between neighboring sites with difi^erent spins, 
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and by independently opening with probability p edges between neighboring sites with 
same spins. Then, for any {u},a), 

Q-^l^r{a)pO(ij)^^ _ p-^a-o(ui)-r(a) pO(u;) (^l _ py(ui) 

P[{^,^)] = —f = —f 

where a is the number of edges of G and r{a) the number of edges between sites with 
different spins. 

Note that the two previous measures are in fact defined on the same set of compatible 
pairs of configurations: if a has been obtained from uj, then uj can be obtained from a 
via the second procedure described above, and the same is true in the reverse direction 
for ui and a. Therefore, P = P and the marginals of P are the FK percolation with 
parameters (p, 2) and the Ising model at inverse-temperature /3, which is the claim. □ 

The coupling gives a randomized procedure to obtain a spin-Ising configuration from 
a FK-Ising configuration (it suffices to assign random spins). The proof of Proposi- 
tion |3.8| provides a randomized procedure to obtain a FK-Ising configuration from a 
spin-Ising configuration. 

If one considers wired boundary conditions for the FK percolation, the Edwards- 
Sokal coupling provides us with an Ising configuration with + boundary conditions (or 
-, the two cases being symmetric). We do not enter into details, since the generalization 
is straightforward. 

An important consequence of the Edwards-Sokal coupling is the relation between 
Ising correlations and FK connectivity properties. Indeed, two sites which are connected 
in the FK percolation configuration must have the same spin, while sites which are not 
have independent spins. This implies: 

Corollary 3.9. For p e (0, 1), G a finite graph and 1^ = -\ lii(l - p), we obtain 

M^j.cK] = (l^l,2,Gi^ ^ cIG). 
In particular, f3c = ~|lii[l -pc(2)]. 

Proof We leave the proof as an exercise. □ 

The uniqueness of Ising infinite- volume measures was discussed in the previous sec- 
tion. The same question can be asked in the case of the FK-Ising model. First, it can be 
proved that (ftp 2 ^^nd (j)^ 2 are extremal among all infinite- volume measures. Therefore, 
it is sufficient to prove that (/>p 2 = '/'p 2 to prove uniqueness. Second, the absence of 
an infinite cluster for (pp 2 can be shown to imply the uniqueness of the infinite- volume 
measure. Using the equality Pc = Psd, the measure is necessarily unique whenever p < p^d 
since (p^ 2 has no infinite cluster. Planar duality shows that the only value of p for which 
uniqueness could eventually fail is the (critical) self-dual point \/2/(l + \/2). It turns 
out that even for this value, there exists a unique infinite volume measure. Since this 
fact will play a role in the proof of conformal invariance, we now sketch an elementary 
proof due to W. Werner (the complete proof can be found in |Wer09] ). 

Proposition 3.10. There exists a unique infinite-volume FK-Ising measure with param- 
eter Pc = \/2/(l + \/2) and there is almost surely no infinite cluster under this measure. 
Correspondingly, there exists a unique infinite-volume spin Ising measure at /3c. 
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Proof As described above, it is sufficient to prove that 4''i,^, 2 = 't']},^ 2- First note that 



there is no infinite cluster for 4>^^^ 2 thanks to Exercise 



3.6 



Via the Edwards-Sokal 



couphng, the infinite- volume Ising measure with free boundary conditions, denoted by 
/j,^ , can be constructed by coloring clusters of the measure cj)^^^ ^. Since there is no 
infinite cluster, this measure is obviously symmetric by global exchange of +/-. In 
particular, the argument of Exercise 3.6 can be applied to prove that there are neither 
+ nor - infinite clusters. Therefore, fixing a box, there exists a + star-connected circuit 
surrounding the box with probability one (two vertices x and y are said to be star- 
connected if y is one of the eight closest neighbors to x). 

One can then argue that the configuration inside the box stochastically dominates the 
Ising configuration for the infinite-volume measure with + boundary conditions (roughly 
speaking, the circuit of spin + behaves like + boundary conditions). We deduce that 
/i^ restricted to the box (in fact to any box) stochastically dominates . This implies 
that /i^ ^ • Since the other inequality is obvious, /i^ and are equal. 

Via Edwards-Sokal's coupling again, 2 = 0^^^ 2 ^"^^ there is no infinite cluster at 
criticality. Moreover, pTp^ = = ^^id there is a unique infinite- volume Ising measure 
at criticality. □ 



Remark 3.11. More generally, the FK percolation with integer parameter q > 2 can 
be coupled with Potts models. Many properties of Potts models are derived using FK 
percolation, since we have the FKG inequality at our disposal, while there is no equivalent 
of the spin-Ising FKG inequality for Potts models. 

3.3 Loop representation of the FK-Ising model and fermionic 
observable 

Let (f2,a,6) be a simply connected domain with two marked points on the boundary. 
Let ris be an approximation of il, and let dab and dba denote the counterclockwise arcs 
in the boundary d^ls joining a to 6 (resp. 6 to a). We consider a FK-Ising measure 
with wired boundary conditions on dba ^ all the edges are pairwise connected - and 
free boundary conditions on the arc dab- These boundary conditions are called the 
Dobrushin boundary conditions. We denote by 'Pq^ p the associated FK-Ising measure 
with parameter p. 

The dual boundary arc d^^ is the set of sites of fig adjacent to dba while the dual 
boundary arc d*ij is the set of sites of L| xfij adjacent to dab, see Fig. [6] A FK-Dobrushin 
domain (flg,as,bs) is given by 

• a medial graph fl^ defined as the set of medial vertices associated to edges of fls 
and to dual edges of 9*^, 

• medial sites as,bs e ftg between arcs dba an 9*f^, see Fig. |6] again, 

with the additional condition that 65 is the southeast corner of a black face belonging 
to the domain. 



Remark 3.12. Note that the definition 0/ fi^ is not the same as in Section 1.2.1 since 
we added medial vertices associated to dual edges of d*f^. We chose this definition to 
make sites of the dual and the primal lattices play symmetric roles. The condition that 
bs is the south corner of a black face belonging to the domain is a technical condition. 

Let {Slg,as,bs) be a FK-Dobrushin domain. For any FK-Ising configuration with 
Dobrushin boundary conditions on fi^, we construct a loop configuration on fig as 
follows: The interfaces between the primal clusters and the dual clusters (i.e clusters in 
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Figure 6: A domain fig with Dobrushin boundary conditions: the vertices of the primal 
graph are black, the vertices of the dual graph fij are white, and between them lies the 
medial graph fig. The arcs dba and 9*;, are the two outermost arcs. Moreover, arcs d^^ 
and dab are the arcs bordering dba and 9*^ from the inside. The arcs dab and dba (resp. 
d*i^ and d^^) are drawn in solid lines (resp. dashed lines) 
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Figure 7: A FK percolation configuration in the Dobrushin domain {Cls,as,bs), together 
with the corresponding interfaces on the medial lattice: the loops are grey, and the 
exploration path 7 from as to bg is black. Note that the exploration path is the interface 
between the open cluster connected to the wired arc and the dual-open cluster connected 
to the white faces of the free arc. 
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the dual model) form a family of loops together with a path from ag to bg . The loops are 
drawn as shown in Figure |7] following the edges of the medial lattice. The orientation 
of the medial lattice naturally gives an orientation to the loops, so that we are working 
with a model of oriented loops on the medial lattice. 

The curve from as to bg is called the exploration path and denoted by 7 = 7(1^). It 
is the interface between the open cluster connected to dba and the dual-open cluster 
connected to 9*^. As in the Ising model case, one can study its scaling limit when the 
mesh size goes to 0: 

Theorem 3.13 (Conformal invariance of the FK-Ising model, [KSIOJ). Let fl be a 
simply connected domain with two marked points a, b on the boundary. Let 75 be the 
interface of the critical FK-Ising with Dobrushin boundary conditions on {D,s,as,bs). 
Then the law of 75 converges weakly, when S ^ 0, to the chordal Schramm- Loewner 
Evolution with k = 16/3. 

As in the Ising model case, the proof of this theorem also involves a discrete observ- 
able, which converges to a conformally invariant object. We define it now. 

Definition 3.14. The edge FK fermionic observable is defined on edges of by 

J^n»,a„.„p(e) = E-/;[e^-'^^(^'''^)le.,], (3.5) 

where W-y{e,bs) denotes the winding between the center of e and bg. 

The vertex FK fermionic observable is defined on vertices of fig \ dfl^ by 

Fn-.ae^bsA^) = lll ^a^,a,,6,,p(e) (3.6) 

where the sum is over the four medial edges having v as an endpoint. 

When we consider the observable at criticality (which will be almost always the case) , 
we drop the dependence on p in the notation. More generally, if (J7,a,5) is fixed, we 
simply denote the observable on {flg,as,bs,Psd) by Fg. 

The quantity ^^(e) is a complexified version of the probability that e belongs to 
the exploration path. The complex weight makes the link between Fg and probabilistic 
properties less explicit. Nevertheless, the vertex fermionic observable Fg converges when 
6 goes to 0: 

Theorem 3.15. jSmilOa^ Let {n,a,b) be a simply connected domain with two marked 
points on the boundary. Let Fg be the verteix. fermionic observable in {Clg,as,bg). Then, 
we have 

h^Fgi-) - V<H^ whenS^O (3.7) 
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uniformly on any compact subset of Q, where (j) is any conformal map from fl to the 
strip M X (0, 1) mapping a to -00 and b to 00. 

As in the case of the spin Ising model, this statement is the heart of the proof of 
conformal invariance. Yet, the observable itself can be helpful for the understanding of 
other properties of the FK-Ising model. For instance, it enables us to prove a statement 
equivalent to the celebrated Russo-Seymour- Welsh Theorem for percolation. This result 
will be central for the proof of compactness of exploration paths (an important step in 



the proof of Theorems 2.10 and 3.13) 



Theorem 3.16 (RSW-type crossing bounds, jDCHNlO] ). There exists a constant c> 
such that for any rectangle R of size An x n, one has 

'Ppsd 2 flC^^^'^^ exists an open path from left to right) > c. (3-8) 
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Before ending this section, we present a simple yet crucial result: we show that it 
is possible to compute rather explicitly the distribution of the loop representation. In 
particular, at criticality, the weight of a loop configuration depends only on the number 
of loops. 

Proposition 3.17. Let p e (0, 1) and let {ilg,as,bs) be a FK Dobrushin domain, then 
for any configuration io, 

= I x°(")x/2'^"' (3.9) 

where x = p/[\/2(l-p)], ^(a;) is the number of loops in the loop configuration associated 
to uj, o{uj) is the number of open edges, and Z is the normalization constant. 



Proof Recall that 



= ^b/(i-p)]°(")2'^("). 



3.3 



the dual of (^^^.'''''i can be proved to be c^^''"-' 



Using arguments similar to Proposition 

(in this sense, Dobrushin boundary conditions are self-dual). With a;* being the dual 
configuration of w, we find 



lbs. as 



o{u>*) I— k{u]*) 



o(w) 



\/ZZ*\ il-p)p* 



p{l-p*) /—TT, — o(uj*) + o{uj) ^k{u)+k{uj*) 

Vp /(l-p ) v2 



2Vp*/(1-P*) 



o{u)) + o{uj* ) 



where the definition of p* was used to prove that = x^. Note that i{uj) 

k{uj) + k{uj*) - 1 and 



2y'p*/{l-p*) 



o(u)+o(lu* ) 



does not depend on the configuration (the sum o(w) + o{uj*) being equal to the total 
number of edges). Altogether, this implies the claim. □ 



4 Discrete complex analysis on graphs 

Complex analysis is the study of harmonic and holomorphic functions in complex do- 
mains. In this section, we shall discuss how to discretize harmonic and holomorphic 
functions, and what are the properties of these discretizations. 

There are many ways to introduce discrete structures on graphs which can be devel- 
oped in parallel to the usual complex analysis. We need to consider scaling limits (as 
the mesh of the lattice tends to zero) , so we want to deal with discrete structures which 
converge to the continuous complex analysis as finer and finer graphs are taken. 
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4.1 Preharmonic functions 

4.1.1 Definition and connection with random walks. 

Introduce the (non-normalized) discretization of the Laplacian operator A 
dyy) in the case of the square lattice L5. For u e L5 and / : L5 ^ C, define 

^^/(") = J E (/(«)- /("))■ 

The definition extends to rescaled square lattices in a straightforward way (for instance 
to L|). 

Definition 4.1. A function h •■ fls ^ C is preharmonic (resp. pre-superharmonic, 
pre-subharmonic j if Ash{x) = (resp. < 0, >0) for every x e fls- 

One fundamental tool in the study of preharmonic functions is the classical relation 
between preharmonic functions and simple random walks: 

Let (Xn) be a simple random walk killed at the first time it exits Vis; then h is 
preharmonic on if and only if {h{Xn)) is a martingale. 

Using this fact, one can prove that harmonic functions are determined by their value 
on dils; that they satisfy Harnack's principle, etc. We refer to |Law91| for a deeper 
study on preharmonic functions and their link to random walks. Also note that the set 
of preharmonic functions is a complex vector space. As in the continuum, it is easy to 
see that preharmonic functions satisfy the maximum and minimum principles. 



4.1.2 Derivative estimates and compactness criteria. 

For general functions, a control on the gradient provides regularity estimates on the 
function itself. It is a well-known fact that harmonic functions satisfy the reverse prop- 
erty: controlling the function allows us to control the gradient. The following lemma 
shows that the same is true for preharmonic functions. 

Proposition 4.2. There exists C > such that, for any preharmonic function h : fls ^ C 
and any two neighboring sites x,y e fig, 

Proof Let x,y e H.^. The preharmonicity of h translates to the fact that h{Xn) is 
a martingale (where Xn is a simple random walk killed at the first time it exits fig). 
Therefore, for x,y two neighboring sites of Qg, we have 

h{x) - h{y) = E[h{Xr) - h{Yr')] (4.2) 

where under E, X and Y are two simple random walks starting respectively at x and 
y, and r, r' are any stopping times. Let 2r = d{x,il.^) > 0, so that U = x + [-r, r]^ is 
included in fig. Fix r and r' to be the hitting times of dUg and consider the following 



coupling of X and Y (one has complete freedom in the choice of the joint law in (4.2 1): 
{X„) is a simple random walk and Yn is constructed as follows, 

• a Xi = y, then Yn = Xn+i for n > 0, 

• if Xi + y, then Yn = a{Xn+i), where a is the orthogonal symmetry with respect to 
the perpendicular bisector £ of [Xi,?/], whenever Xn+i does not reach i. As soon 
as it does, set y„ = Xn+i- 
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It is easy to check that Y is also a simple random walk. Moreover, we have 

\h{x)-h{y)\<E[\h{Xr)-h{Yr')\lx^^Y^,]<2\ Snp \h{z)\] V{Xr+Yr') 

Using the definition of the coupling, the probability on the right is known: it is equal 
to the probability that X does not touch t before exiting the ball and is smaller than 
^(5 (with C a universal constant), since Us is of radius rjd for the graph distance. We 
deduce that 

|/i(a:)-/i(y)|<2( sup \h{z)\\—5 < 2 I sup \h{z)\]— 5 
\zedUs I r I r 

a 



Recall that functions on fls are implicitly extended to fl. 

Proposition 4.3. A family {hs)s>a of preharmonic functions on the graphs fls is pre- 
compact for the uniform topology on compact subsets of Q if one of the following prop- 
erties holds: 

(1) {hs)syQ is uniformly bounded on any compact subset ofVt, 

or 

(2) for any compact subset K of fJ, there exists M = M{K) > such that for any 
S>0, 

5^ ^ \hs{x)\''<M. 

x^Ks 



Proof Let us prove that the proposition holds under the first hypothesis and then that 
the second hypothesis implies the first one. 

We are faced with a family of continuous maps kg : D, ^ C and we aim to apply 
the Arzela-Ascoli theorem. It is sufficient to prove that the functions hs are uniformly 
Lipschitz on any compact subset since they are uniformly bounded on any compact 
subset of f2. Let K he a, compact subset of il. Proposition 4.2 shows that \hs{x)-hs{y)\ < 
Ck5 for any two neighbors x, y e Kg, where 

_ SUPa>0 ^^Vxen:d{x.K)<rl2 \hs{x)\ 

implying that \hs{x) - hs{y)\ < 2Ck\x - y\ for any x,y e Kg (not necessarily neighbors). 
The Arzela-Ascoli theorem concludes the proof. 

Now assume that the second hypothesis holds, and let us prove that {hs)s>o is 
bounded on any compact subset of il. Take K cQ, compact, let 2r = d{K^Vl'^) > and 
consider x e Kg. Using the second hypothesis, there exists k ■■= k{x) such that ^ < k < ^ 
and 

S E \''s{y)\' < 2M/r, (4.3) 

yedUks 

where Uks = x + [-Skjdk^^ is the box of size k (for the graph distance) around x and 
M = M{y + [-r, r]^). Exercise 4.4 implies 



hs{x) = E hs{y)Hu^,{x,y) (4.4) 
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for every x € Usk- Using the Cauchy-Schwarz inequality, we find 

\yeaUks I 

where C is a uniform constant. The last inequality used Exercise |4.5| to affirm that 
Huks {^^ v) ^ C'cJ for some C = C{r) > 0. □ 



Exercise 4.4. The discrete harmonic measure HQ^{-,y) of y e dQs is the unique har- 
monic function on \dVLs vanishing on the boundary Oils, except at y, where it equals 
1. Equivalently, Hfi^(x,y) is the probability that a simple random walk starting from x 
exits D.S \ dils through y. Show that for any harmonic function /i : fl^ -> C, 

h = E Hy)HnA;y)- 

Exercise 4.5. Prove that there exists C > such that Hq^ (0, y) < C5 for every 5 > 
and y € dQs, where Q = [-1, 1]^. 



4.1.3 Discrete Dirichlet problem and convergence in the scaling limit. 

Preharmonic functions on square lattices of smaller and smaller mesh size were studied 
in a number of papers in the early twentieth century (see e.g. |PW23I rBou26i lLus26| l. 
culminating in the seminal work of Courant, Friedrichs and Lewy. It was shown in 
[CFL28| that solutions to the Dirichlet problem for a discretization of an elliptic operator 
converge to the solution of the analogous continuous problem as the mesh of the lattice 
tends to zero. A first interesting fact is that the limit of preharmonic functions is indeed 
harmonic. 

Proposition 4.6. Any limit of a sequence of preharmonic functions on fig converging 
uniformly on any compact subset of is harmonic in Q. 



Proof Let (hs) be a se quence of preharmonic functions on ils converging to h. Via 



Propositions 



4.2 



and 



4.3 



{■T[hs{- + S) - hs])s>o is precompact. Since dxh is the only 



possible sub-sequential limit of the sequence, {^^[iT-si' + S) - hs])s>o converges (indeed 
its discrete primitive converges to h) . Similarly, one can prove convergence of discrete 
derivatives of any order. In particular, = -^Aghg converges to ^[dxxh + dyyh]. There- 
fore, h is harmonic. □ 



In particular, preharmonic functions with a given boundary value problem converge 
in the scaling limit to a harmonic function with the same boundary value problem in a 
rather strong sense, including convergence of all partial derivatives. The finest result of 
convergence of discrete Dirichlet problems to the continuous ones will not be necessary 
in our setting and we state the minimal required result: 

Theorem 4.7. Let CI be a simply connected domain with two marked points a and 
b on the boundary, and f a bounded continuous function on the boundary of f2. Let 
fs '■ dfls ->■ C be a sequence of uniformly bounded functions converging uniformly away 
from a and b to f. Let hs be the unique preharmonic map on fls such that {hs)\Qn^ = fs- 
Then 

hs — > h when S ^ 
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uniformly on compact subsets of Q, where h is the unique harmonic function on Q, 
continuous on fl, satisfying h^g^ = f . 

Proof Since {fs)5>a is uniformly bounded by some constant M, the minimum and 
maximum principles imp ly th at {hs)s>o is bounded by M. Therefore, the family (hs) is 
precompact (Proposition 4.3 1. L et h be a sub-sequential limit. Necessarily, /i is harmonic 
inside the domain (Proposition 4.6 1 and bounded. To prove that h = h, it suffices to 
show that h can be continuously extended to the boundary by /. 

Let X € d^l \ {a, b} and e > 0. There exists R> such that for S small enough, 

\fs{x') - f5{x)\ < £ for every x' edVln Q{x, R), 

where Q{x, R) = x + [-R, R]^. For r < R and y e Q{x, r), we have 

\hs{y) - fs{x)\ = Ey[fs{Xr) - fs{x)] 

for X a random walk starting at y, and r its hitting time of the boundary. Decomposing 
between walks exiting the domain inside Q{x,R) and others, we find 

\hs{y) - f5{x)\ < e + 2MFy[XriQ{x,R)] 

Exercise 4.8 guarantees that Py[XT- ^ Q{x, R)] < {rjR)'^ for some independent constant 
a > 0. Taking r = R{e/2My/" and letting S go to 0, we obtain \h{y) - f{x) \ < 2e for 
every y e Q{x, r). □ 



Exercise 4.8. Show that there exists a > such that for any 1 » r > 6 > and any 
curve 7 inside D := {z : \z\ < 1} from C = {z : \z\ = 1} to {z : \z\ = r}, the probability for a 
random walk on starting at to exit (D\j)s through C is smaller than r" . To prove 
this, one can show that in any annulus {z : a; < |z| < 2x}, the random walk trajectory has 
a uniformly positive probability to close a loop around the origin. 

4.1.4 Discrete Green functions 

This paragraph concludes the section by mentioning the important example of discrete 
Green functions. For j/ e 57^ \ dfls, let Gfi^{-,y) be the discrete Green function in the 
domain Qg with singularity at y, i.e. the unique function on fls such that 

• its Laplacian on rig \ dfls equals except at y, where it equals 1, 

• (•;?/) vanishes on the boundary dflg. 

The quantity -Gq^{x, y) is the number of visits at a; of a random walk started at y and 
stopped at the first time it reaches the boundary. Equivalently, it is also the number of 
visits at y of a random walk started at x stopped at the first time it reaches the boundary. 
Green functions are very convenient, in particular because of the Riesz representation 
formula for (not necessarily harmonic) functions: 

Proposition 4.9 (Riesz representation formula). Let f •■ fls ^ C be a function vanishing 
on dQg. We have 

f = Z A,/(y)Go,(-,y). 
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Proof Note that /- ZyeOj '^sf{y)Gns {'i v) is harmonic and vanishes on the boundary. 
Hence, it equals everywhere. □ 



Finally, a regularity estimate on discrete Green functions will be needed. This 
proposition is slightly technical. In the following, aQs = [-a,a]^ n L5 and Vxfix) = 
{f{x + S)-f{x),f{x + zS)-f{x)). 

Proposition 4.10. There exists C > such that for any S > and y e 9Qs, 
E |VxG9Q,(x,y)| < CS GgQ,{x,y). 

Proof In the proof, Ci,...,Cq denote universal constants. First assume y e 9Qs \ SQs- 
Using random walks, one can easily show that there exists Ci > such that 

■^G9Q,{x,y) < G9Q,{x',y) < CiGgg, (x, ?/) 

for every x, x' e 2Qs (this is a special application of Harnack's principle). Using Propo- 
sition |42j we deduce 

Y, WxGgQ,{x,y)\< G2S maxGgQ,{x,y) <GiC2S Y GgQ,{x,y) 

xeQs xeQs xi2Qs ^^q^ 

which is the claim for y e 9Qs \ iQg. 

Assume now that y e iQg. Using the fact that Ggg^ (a;, y) is the number of visits of 
X for a random walk starting at y (and stopped on the boundary), we find 

Y G9Qsix,y)>C,/S'. 

Therefore, it suffices to prove HxnQi, I VGgg^ (x, y)\ < C4/S. Let Gl^ be the Green function 
in the whole plane, i.e. the function with Laplacian equal to Sx^y, normalized so that 
Ghs (z/j y) = 0, and with sublinear growth. This function has been widely studied, it was 
proved in |MW40| that 



GL,(a;,2/) = 




Now, Gl^ ('!?/)" Ggg^, (•, J/) - ;^ In is harmonic and has bounded boundary conditions 
on d9Qs. Therefore, Proposition |4.2| implies 

Y |v.(GL,(x,y)-G9Q,(x,y))| < C^S-l/S^ = Ge/S. 

Moreover, the asymptotic of Gl^ (■,?/) leads to 

Y |v^GL,(a;,y)| < Gt/^. 

xnQs 

Summing the two inequalities, the result follows readily. □ 
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4.2 Preholomorphic functions 



4.2.1 Historical introduction 

Preholomorphic functions appeared implicitly in KirchhofF's work |Kir47| . in which a 
graph is modeled as an electric network. Assume every edge of the graph is a unit 
resistor and for u ~ v, let F{uv) be the current from u to v. The first and the second 
KirchhofF's laws of electricity can be restated: 

• the sum of currents flowing from a vertex is zero: 

E = 0' (4-5) 



• the sum of the currents around any oriented closed contour 7 is zero: 

E F{uv) = 0. (4.6) 



Different resistances amount to putting weights into (4.51 and (4.6 1. The second law 
is equivalent to saying that F is given by the gradient of a potential function H, and 
the first equivalent to H being preharmonic. 

Besides the original work of Kirchhoff, the first notable application of preholomorphic 
functions is perhaps the famous article |BSST40] of Brooks, Smith, Stone and Tutte, 
where preholomorphic functions were used to construct tilings of rectangles by squares. 

Preholomorphic functions distinctively appeared for the first time in the papers 
|Isa41i IIsa52| of Isaacs, where he proposed two definitions (and called such functions 
mono-diffric) . Both definitions ask for a discrete version of the Cauchy-Riemann equa- 
tions diaF = idaF or equivalently that the z-derivative is 0. In the first definition, the 
equation that the function must satisfy is 

i[f{E)-f{S)] = f{W)-f{S) 

while in the second, it is 

^[f{E)-f{W)] = f{N)-f{S), 

where N, E, S and W are the four corners of a face. A few papers of his and other 
mathematicians followed, studying the first definition, which is asymmetric on the square 
lattice. The second (symmetric) definition was reintroduced by Ferrand, who also dis- 
cussed the passage to the scaling limit and gave new proofs of Riemann uniformization 
and the Courant-Friedrichs-Lewy theorems per44. .LF55j . This was followed by exten- 
sive studies of Duflin and others, starting with |Duf56| . 



4.3 Isaacs's definition of preholomorphic functions 

We will be working with Isaacs's second deflnition (although the theories based on both 
deflnitions are almost the same). The deflnition involves the following discretization of 
the d = dx + idy operator. For a complex valued function / on hg (or on a flnite subgraph 
of it) , and x e , deflne 

dsfix) = l[fiE)-fiW)] + '-[fiN)-fiS)] 

where N, E, S and W denote the four vertices adjacent to the dual vertex x indexed in 
the obvious way. 
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Remark 4.11. When defining derivation, one uses duality between a graph and its dual. 
Quantities related to the derivative of a function on G are defined on the dual graph G* . 
Similarly, notions related to the second derivative are defined on the graph G again, 
whereas a primitive would he defined on G* . 

Definition 4.12. A function f : ilg ^ C is called preholomorphic if dsf{x) = for 
every x e fig. For x e il*g, dsf{x) = is called the discrete Cauchy-Riemann equation at 

X. 

The theory of preholomorphic functions starts much Hke the usual complex analysis. 
Sums of preholomorphic functions are also preholomorphic, discrete contour integrals 
vanish, primitive (in a simply-connected domain) and derivative are well-defined and are 
preholomorphic functions on the dual square lattice, etc. In particular, the (discrete) 
gradient of a preharmonic function is preholomorphic (this property has been proposed 
suitable generalization in higher dimensions). 

Exercise 4.13. Prove that the restriction of a continuous holomorphic function to hs 
satisfies discrete Cauchy-Riemann equations up to 0{5^). 

Exercise 4.14. Prove that any preholomorphic function is preharmonic for a slightly 
modified Laplacian (the average over edges at distance \/2S minus the value at the point). 
Prove that the (discrete) gradient of a preharmonic function is preholomorphic (this 
property has been proposed as a suitable generalization in higher dimensions). Prove 
that the limit of preholomorphic functions is holomorphic. 

Exercise 4.15. Prove that the integral of a preholomorphic function along a discrete 
contour vanishes. Prove that the primitive and the differential of preholomorphic func- 
tions are preholomorphic. 

Exercise 4.16. Prove that -^j^ds and -^^s converge (when 5^0) to d, B and A in 
the sense of distributions. 

4.4 s-holomorphic functions 

As explained in the previous sections, the theory of preholomorphic functions starts like 
the continuum theory. Unfortunately, problems arrive quickly. For instance, the square 
of a preholomorphic function is no longer preholomorphic in general. This makes the 
theory of preholomorphic functions significantly harder than the usual complex analysis, 
since one cannot transpose proofs from continuum to discrete in a straightforward way. 
In order to partially overcome this difficulty, we introduce s-holomorphic functions (for 
spm-holomorphic), a notion that will be central in the study of the spin and FK fermionic 
observables. 

4.4.1 Definition of s-holomorphic functions 

To any edge of the medial lattice e, we associate a line £(e) passing through the origin 
and \/l (the choice of the square root is not important, and recall that e being oriented, 
it can be thought of as a complex number) . The different lines associated with medial 
edges on are M, e^Z-^M, iR and e^'^^/^M, see Fig. [s] 

Definition 4.17. A function f:Clg^C is s-holomorphic if for any edge e ofQg, we 
have 

Pi(e)[m] = Pi(e)[f{y)] 

where x,y are the endpoints of e and Pp is the orthogonal projection on I. 
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Figure 8: Lines £{e) for medial edges around a white face. 

The definition of s-holomorphicity is not rotationally invariant. Nevertheless, / is 
s-holoniorphic if and only if e*'^^*/(z-) (resp. «/(-•)) is s-holomorphic. 

Proposition 4.18. Any s-holomorphic function / : il^ C is preholomorphic on fig. 

Proof Let / : fi^ ^ C be a s-holomorphic fimction. Let w be a vertex of hg u (this 
is the vertex set of the dual of the medial lattice). Assume that v e ilg, the other case 
is similar. We aim to show that dsf{v) = 0. Let NW, NE, SE and SW be the four 
vertices around v as illustrated in Fig. |8] Next, let us write relations provided by the 
s-holomorphicity, for instance 

Pk[/(W)] = P^[f{NE)l 

Expressed in terms of / and its complex conjugate / only, we obtain 

f{NW)+f{NW) = f{NE)+J{NE). 

Doing the same with the other edges, we find 

f{NE)+ iJ{NE ) = fiSE)+i JiSE) 
fiSE)-f iSE) = fiSW)-f iSW) 
f{SW)-if{SW) = f{NW)-if(NW) 

Multiplying the second identity by -i, the third by -1, the fourth by i, and then summing 
the four identities, we obtain 

0={l-i) [f{NW) - f{SE) + ifiSW) - tf{NE)] = 2(1 - t)dsf{v) 

which is exactly the discrete Cauchy-Riemann equation in the medial lattice. □ 

4.4.2 Discrete primitive of 

One might wonder why s-holomorphicity is an interesting concept, since it is more 
restrictive than preholomorphicity. The answer comes from the fact that a relevant 
discretization of ilm(/^/^) can be defined for s-holomorphic functions /. 

Theorem 4.19. Let f : Qg ^ C be an s-holomorphic function on the discrete simply 
connected domain fig, and bo ^ fls- Then, there exists a unique function H •■ nguflg -> C 
such that 

H(b{)) = 1 and 
H{b)-H{w) = <5 |F,(e)[/(a;)]r (= <5 |F,(e)[/(y)]r ) 
for every edge e = [xy] of fig bordered by a black face b e fig and a white face w e fig. 
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Figure 9: Arrows corresponding to contributions to 2AH' . Note that arrows from black 
to white contribute negatively, those from white to black positively. 



An elementary computation shows that for two neighboring sites 61,62 e i^s, with v 
being the medial vertex at the center of [6162], 

H{h)-H{b2) = hm[f{vf-{h-b2)], 

the same relation holding for sites of J . This legitimizes the fact that iJ is a discrete 
analogue of |lm (/^ /^). 

Proof The uniqueness of H is straightforward since Qg is simply connected. To obtain 
the existence, construct the value at some point by summing increments along an arbi- 
trary path from 60 to this point. The only thing to check is that the value obtained does 
not depend on the path chosen to define it. Equivalently, we must check the second 
Kirchhoff's law. Since the domain is simply connected, it is sufficient to check it for 
elementary square contours around each medial vertex v (these are the simplest closed 
contours). Therefore, we need to prove that 

\Pi(n)[f{v)f - \Piie)[f{v)f + \Piis)[f{v)f " \Pii^u)[f {v)f = 0, (4.7) 

where n, e, s and w are the four medial edges with endpoint v, indexed in the obvious 



way. Note that £{n) and £(s) (resp. £{e) and £{w)) are orthogonal. Hence, (4.7 1 follows 
from 

\Piin)[fiv)f + \P,is)[fiv)f = \fiv)\' = \Piie)[fiv)f + \Pii^^[fiv)f. (4.8) 



Even if the primitive of / is preholomorphic and thus preharmonic, this is not the case 
for H in generaj^ Nonetheless, H satisfies subharmonic and superharmonic properties. 
Denote by H' and H° the restrictions of _ff : 17^ u fij C to fig (black faces) and 
(white faces). 

Proposition 4.20. If f : ^ C is s-holomorphic, then H' and H° are respectively 
subharmonic and superharmonic. 

^ H is roughly (the imaginary part of) the primitive of the square of /. 
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Figure 10: The black graph is tlie isoradial graph. Grey vertices are the vertices on the 
dual graph. There exists a radius r > such that all faces can be put into an incircle 
of radius r. Dual vertices have been drawn in such a way that they are the centers of 
these circles. 

Proof Let B be a vertex of ^ d^s- We aim to show that the sum of increments of 
H' between B and its four neighbors is positive. In other words, we need to prove that 
the sum of increments along the sixteen arrows drawn in Fig. |9]is positive. Let a, 6, c 
and d be the four values of \/Jff(e)[/(2/)] for every vertex y e 17^ around B and any 
edge e = [yz\ bordering B (there are only four different values thanks to the definition 
of s-holomorphicity) . An easy computation shows that the eight interior increments are 
thus -a^, -6^, -(?, -(P (each appearing twice). Using the s-holomorphicity of / on 
vertices of fig around B, we can compute the eight exterior increments in terms of a, b, 
cand d: we obtain (a^/2-6)^ {bV2-af, {bV2-cf, {c^/2-bY, {cV^-df, (dV^-c)^, 
(d\/2 + a)2, (a\/2 + d)2. Hence, the sum S of increments equals 



Remark 4.21. A subharmonic function in a domain is smaller than the harmonic 
function with the same boundary conditions. Therefore, H' is smaller than the harmonic 
function solving the same boundary value problem while H° is bigger than the harmonic 
function solving the same boundary value problem. Moreover, H'(b) is larger than 
H°(w) for two neighboring faces. Hence, if H' and H° are close to each other on the 
boundary, then they are sandwiched between two harmonic functions with roughly the 
same boundary conditions. In this case, they are almost harmonic. This fact will be 
central in the proof of conformal invariance. 

4.5 Isoradial graphs and circle packings 

Duffin [Duf68| extended the definition of preholomorphic functions to isoradial graphs. 
Isoradial graphs are planar graphs that can be embedded in such a way that there exists 



S = A{a^ + + c^ + d^)-4\/2{ab + bc + cd-da) 
= 4|e-"/*a-5 + e*3"/'*c-id|' > 0. 



(4.9) 
(4.10) 



The proof for H° follows along the same lines. 



□ 
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r > so that each face has a circumcircle of same radius r > 0, see Fig. [TO] When the 
embedding satisfies this property, it is said to be an isoradial embedding. We would 
like to point out that isoradial graphs form a rather large family of graphs. While not 
every topological quadrangulation (graph all of whose faces are quadrangles) admits a 
isoradial embedding, Kenyon and Schlenker |KS05| gave a simple necessary and sufficient 
topological condition for its existence. It seems that the first appearance of a related 
family of graphs in the probabilistic context was in the work of Baxter |Bax89| , where the 
eight- vertex model and the Ising model were considered on Z-invariant graphs, arising 
from planar line arrangements. These graphs are topologically the same as the isoradial 
ones, and though they are embedded differently into the plane, by |KS05| they always 
admit isoradial embeddings. In |Bax89| . Baxter was not considering scaling limits, and 
so the actual choice of embedding was immaterial for his results. However, weights in his 
models would suggest an isoradial embedding, and the Ising model was so considered by 
Mercat |Mer01| . Boutilier and de Tihere |BdTlltlBdTT0l . Chelkak and Smirnov [CSOS] 
(see the last section for more details). Additionally, the dimer and the uniform spanning 
tree models on such graphs also have nice properties, see e.g. |Ken02| . Today, isoradial 
graphs seem to be the largest family of graphs for which certain lattice models, including 
the Ising model, have nice integrability properties (for instance, the star-triangle relation 
works nicely). A second reason to study isoradial graphs is that it is perhaps the largest 
family of graphs for which the Cauchy-Riemann operator admits a nice discretization. 
In particular, restrictions of holomorphic functions to such graphs are preholomorphic 
to higher orders. The fact that isoradial graphs are natural graphs both for discrete 
analysis and statistical physics sheds yet another light on the connection between the 
two domains. 

In |Thu86| . Thurston proposed circle packings as another discretization of complex 
analysis. Some beautiful applications were found, including yet another proof of the 
Riemann uniformization theorem by Rodin and Sullivan [RS87| . More interestingly, cir- 
cle packings were used by He and Schramm |HS93| in the best result so far on the Koebe 
uniformization conjecture, stating that any domain can be conformally uniformized to 
a domain bounded by circles and points. In particular, they established the conjecture 
for domains with countably many boundary components. More about circle packings 
can be learned from Stephenson's book |Ste05| . Note that unlike the discretizations 
discussed above, the circle packings lead to non-linear versions of the Cauchy-Riemann 
equations, see e.g. the discussion in |BMS05| . 



5 Convergence of fermionic observables 



In this section, we prove the convergence of fermionic observables at criticality (Theo- 
2.11 and 3.151. We start with the easier case of the FK-Ising model. We present 



rems 



the complete proof of the convergence, the main tool being the discrete complex analysis 
that we developed in the previous section. We also sketch the proof of the convergence 
for the spin Ising model. 



5.1 Convergence of the FK fermionic observable 

In this section, fix a simply connected domain (Q, a, b) with two points on the boundary. 
For S > 0, always consider a discrete FK Dobrushin domain {ng,as,bs) and the critical 
FK-Ising model with Dobrushin boundary conditions on it. Since the domain is fixed, 
set Fs = Ffi^ ag.bs.psd ^"-"^ ^^'^ -^-^ fermionic observable. 
The proof of convergence is in three steps: 

• First, prove the s-holomorphicity of the observable. 
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• Second, prove the convergence of the function Hg naturally associated to the s- 
holomorphic functions Fs/\/2d. 

• Third, prove that Fs/\/2S converges to \/(^. 

5.1.1 s-holomorphicity of the (vertex) fermionic observable for FK-Ising. 

The next two lemmata deal with the edge fermionic observable. They are the key steps 
of the proof of the s-holomorphicity of the vertex fermionic observable. 

Lemma 5.1. For an edge e e fi^, Fs{e) belongs to £{e). 

Proof The winding at an edge e can only take its value in the set W + 27rZ where W 
is the winding at e of an arbitrary interface passing through e. Therefore, the winding 
weight involved in the definition of ^^(e) is always proportional to e''^^^ with a real 
coefficient, thus ^^(e) is proportional to e*^^^. In any FK Dobrushin domain, bs is the 
southeast corner and the last edge is thus going to the right. Therefore e'^^^ belongs 
to £{e) for any e and so does Fs{e). □ 

Even though the proof is finished, we make a short parenthetical remark: the defi- 
nition of s-holomorphicity is not rotationally invariant, nor is the definition of FK Do- 
brushin domains, since the medial edge pointing to bg has to be oriented southeast. The 
latter condition has been introduced in such a way that this lemma holds true. Even 
though this condition seems arbitrary, it has no infiuence on the convergence result, 
meaning that one could perform a (slightly modified) proof with another orientation. 

Lemma 5.2. Consider a medial vertex v in fig \ dfl"^. We have 

Fs{N) - Fs{S) = i[Fs{E) - Fs{W)] (5.1) 
where N , E, S and W are the adjacent edges indexed in clockwise order. 



Proof Let us assume that v corresponds to a primal edge pointing SE to NW , see 
Fig. [11] The case NE to SW is similar. 

We consider the involution s (on the space of configurations) which switches the 
state (open or closed) of the edge of the primal lattice corresponding to v. Let e be an 
edge of the medial graph and set 



££7 



the contribution of the configuration u to Fs{e). Since s is an involution, the following 
relation holds: ^ 



In order to prove (5.1 1, it suffices to prove the following for any configuration cj: 

Nu. + iV,(^) - 5^ - = i[E^ + -W^- (5.2) 
There are three possibilities: 

Case 1: the exploration path 7(0;) does not go through any of the edges adjacent to v. 
It is easy to see that neither does 7(s(w)). All the terms then vanish and (5.2) trivially 
holds. 

Case 2: 7(aj) goes through two edges around v. Note that it follows the orientation of 
the medial graph, and thus enters v through either W 01 E and leaves through N or S. 
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Figure 11: Two associated configurations, one with one exploration path and a loop, 
one without the loop. One can go from one to the other by switching the state of the 
edge. 



We assume that 7(w) enters through the edge W and leaves through the edge S {i.e. 
that the primal edge corresponding to v is open) . The other cases are treated similarly. 
It is then possible to compute the contributions of all the edges adjacent to w of w and 
s{uj) in terms of W^j. Indeed, 

• The probability of s{uj) is equal to l/\/2 times the probability of ui (due to the 
fact that there is one less open edge of weight 1 - we are at the self-dual point - 



and one less loop of weight \/2, see Proposition 3.17) 



• Windings of the curve can be expressed using the winding of W. For instance, 
the winding of N in the configuration to is equal to the winding of W minus a 7r/2 
turn. 

The contributions are given as: 



configuration 


W 


E 


N 


S 














s{uj) 






e-W4vK„/y2 





Using the identity e^'^l^ - e~^^/^ = i\/2, we deduce ( |5.2[ ) by summing (with the right 
weight) the contributions of all the edges around v. 

Case 3: 7(0;) goes through the four medial edges around v. Then the exploration path 
of s{oj) goes through only two, and the computation is the same as in the second case. 
In conclusion, 



5.2 1 is always satisfied and the claim is proved. 



Recall that the FK fermionic observable is defined on medial edges as well as on 
medial vertices. Convergence of the observable means convergence of the iiertez observ- 
able. The edge observable is just a very convenient tool in the proof. The two previous 
properties of the edge fermionic observable translate into the following result for the 
vertex fermionic observable. 



37 



Proposition 5.3. The vertex fermionic observable Fs is s-holomorphic. 



Proof Let w be a medial vertex and let N, E, S and W be the four medial edges 



around it. Using Lemmata |5.1| and 5.2 one can see that (5.11 can be rewritten (by 
taking the complex conjugate) as: 

Fs{N)+FsiS) = FsiE) + FsiW). 

In particular, from (|3.6|), 



Fs{v) Pde) = Fs{N)+Fs{S) = FsiE)+Fs{W). 

e adjacent 



Using Lemma 5.1 again, Fs{N) and Fs{S) are orthogonal, so that Fs{N) is the pro- 
jection of Fs{v) on £{N) (and similarly for other edges). Therefore, for a medial edge 
e = [xy], Fs{e) is the projection of Fs{x) and Fs{y) with respect to £{e), which proves 
that the vertex fermionic observable is s-holomorphic. □ 



The function Fs/\/2S is preholomorphic for every S > 0. Moreover, Lemma 5.1 
identifies the boundary conditions of Fg / \/25 (its argument is determined) so that this 
function solves a discrete Riemann-Hilbert boundary value problem. These problems 
are significantly harder to handle than the Dirichlet problems. Therefore, it is more 
convenient to work with a discrete analogue of liJi(^f^[Fs{z)/\/26]^dz^, which should 
solve an approximate Dirichlet problem. 

5.1.2 Convergence of {Hs)s>o. 

Let A be the black face (vertex of fls) borderi ng a/; , see Fig. |6] Since the FK fermionic 



observable Fs/\/2S is s-holomorphic. Theorem 4.19 defines a function Hg : flgLiflg 
such that 



Hs{A) = 1 and 



2 



Hs{B)-Hs{W) = \Peie)[Fs{x)]\ = \Piie)[F5{y)]\ 

for the edge e = [xy] of fig bordered by a black face B e ^Ig and a white face W e 
n*g. Note that its restriction H' to ils is subharmonic and its restriction Hg to ilj is 
superharmonic . 

Let us start with two lemmata addressing the question of boundary conditions for 

Hs- 

Lemma 5.4. The function H' is equal to 1 on the arc dta- The function Hg is equal 
to on the arc 5*^,. 



Proof We first prove that i/* is constant on 9ba- Let B and B' be two adjacent 
consecutive sites of dba- They are both adjacent to the same dual vertex W eUX, see 



Fig. 12 Let e (resp. e') be the edge of the medial lattice between W and B (resp. B'). 



We deduce 



Hl{B)-Hl{B') = \Fs{e)\^-\Fs{e't = 



(5.3) 



The second equality is due to |-F)5(e)| = (j)^<>^^ {W ^ d*^^) (see Lemma 



is constant along the arc. Since Hl{A) = 1, the result follows readily. 



7.31. Hence, m 
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Figure 12: Two adjacent sites B and B' on dha together with the notation needed in the 
proof of Lemma |5.4[ 



Similarly, Hg is constant on the arc Moreover, the dual white face A* e d^i^ 
bordering as (see Fig. [61 satisfies 



H^Ml = Hl{A)-\Fs{e) 



1-1 







(5.4) 



where e is the edge separating A and A* , which necessarily belongs to 7. Therefore 



m = on d:,. 



Lemma 5.5. The function H' converges to on the arc dab uniformly away from a 
and b, Hg converges to 1 on the arc d^^ uniformly away from a and b. 



Proof Once again, we prove the result for H'. The same reasoning then holds for Hg. 
Let i? be a site of dab at distance r of dba (and therefore at graph distance r/S of dba in 
ils)- Let W be an adjacent site of B on 9*^. Lemma 5.4 implies Hg{W) = 0. From the 
definition of Hs, we find 



HUB) = Hl{W) + \P,^M{e)]\ = \Pife)[Fs{e)] 



,(e6 7)' 



Note that e e 7 if and only if B is connected to the wired arc di,a- Therefore, 4>'^^''p ^(e ^ 
7) is equal to the probability that there exists an open path from B to dba (the winding 



is deterministic, see Lemma 7.3 for details). Since the boundary conditions on dab are 
free, the comparison between boundary conditions shows that the latter probability 
is smaller than the probability that there exists a path from B to dUs in the box 
Us = {B + [-J^, r]^) nhs with wired boundary conditions. Therefore, 



H'siB) 



Proposition 3.10 implies that the right hand side converges to (there is no infinite 
cluster for (j^p^^ 2) i which gives a uniform bound for B away from a and b. □ 

The two previous lemmata assert that the boundary conditions for H' and Hg are 
roughly on the arc dab and 1 on the arc dba- Moreover, H' and Hg are almost 
harmonic. This should imply that {Hs)5>q converges to the solution of the Dirichlet 
problem, which is the subject of the next proposition. 
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Proposition 5.6. Let {Q,a,b) be a simply connected domain with two points on the 
boundary. Then, {Hs)s>a converges to Im{(j)) uniformly on any compact subsets of D, 
when d goes to 0, where (j) is any conformal map from Q to T = M x (0, 1) sending a to 
-oo and b to oo. 

Before starting, note that Ini((/)) is the solution of the Dirichlet problem on (51, a, 6) 
with boundary conditions 1 on dba and on dab- 



Proof From the definition of H, Hg is subharmonic. let kg be the preharmonic function 



with same boundary conditions as Hg on dQs 



Note that iJ* < h'g. 



Similarly, kg is defined 



to be the preharmonic function with same boundary conditions as Hg on dQg . If K c fl 
is fixed, where K is compact, let bg e Kg and wg e Kg any neighbor of bg, we have 



h°g{wg) < Hiiwg) < H-g{bg) < h-g{bg). 



(5.5) 



Using Lemmata 5.4 and 5.5 boundary conditions for H' (and therefore h') are uniformly 
and 1 on dba away from a and b. Moreover, \h'g\ is bounded 

kg converges to Im((/)) on 
The 



4.7 



converging to on dab 

by 1 everywhere. This is sufficient to apply Theorem 

any compact subset of when 6 goes to 0. The same reasoning applies to kg. 
convergence for H' and Hg follows easily since they are sandwiched between h' and 



5.1.3 Convergence of FK fermionic observables {Fg/\/2S)gyo. 

This section contains the proof of Theorem |3.15| The strategy is straightforward: 
{Fg I \/26) gyQ is proved to be a precompact family for the uniform convergence on com- 
pact subsets of Q.. Then, the possible sub-sequential limits are identified using Hg. 



Proof of Theorem 3.15 First assume that the precompactness of the family {Fg/\/26)g^Q 
has been proved. Let {Fg^/\/25n)n(:N be a convergent subsequence and denote its limit 
by /. Note that / is holomorphic as it is a limit of preholomorphic functions. For two 
points x,y eil, we have: 

Hg^{y)-Hg„{x) = ^Fliz)dz) 

(for simplicity, also denote the closest points of x,y in flg^^ by x,y). On the one hand, 
the convergence of {Fg^l \/25n)n<i¥i being uniform on any compact subset of Jl, the right 
hand side converges to Im (/J f{z)^dz). On the other hand, the left-hand side converges 
to Im((/)(?/) - (j){x)). Since both quantities are holomorphic functions of y, there exists 
C € K such that (j){y) - 4>{x) = C + f{z)^dz for every x,y eV,. Therefore / equals 
Since this is true for any converging subsequence, the result follows. 

Therefore, the proof boils down to the precompactness of (Fg f \/26) g>Q . We will use 
the second criterion in Proposition |4.3| Note that it is sufficient to prove this result for 
squares Q c n such that a bigger square 9Q (with same center) is contained in n. 

Fix S > 0. When jumping diagonally over a medial vertex v, the function Hg changes 
by Ke{Fg (v)) or lm{Fg{v)) depending on the direction, so that 



Z \Fg{v)/V2Sf = S \VH'g{x)\ + S Z \yHg{x)\ (5.6) 



<;2 



where VH'g{x) = {H'g{x + 5)- H'g{x),H'g{x + iS) - Hg{x)), and \/Hg is defined similarly 
for Hg. It follows that it is enough to prove uniform boundedness of the right hand side 
in (5.6). We only treat the sum involving H'. The other sum can be handled similarly. 
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Write H' = Ss + Rs where Ss is a harmonic function with the same boundary con- 
ditions on d9Qs as H'. Note that Rg < is automatically subharmonic. In order to 
prove that the sum of \VH'\ on Qs is bounded by C/5, we deal separately with \vSs\ 
and |Vi?5|. First, 



E \ySs{x)\ < ^-C^slsnp |^,(x)|]<^( sup \H-s{x)\] 

xeQ, \xeaQs / \xe9Qs / 



- 6 ' 



where in the first inequality we used Proposition |4.2| and the maximum principle for 
Sg, and in the second the fact that 5*^ and H' share the same boundary conditions on 
9Qs- The last inequality comes from the fact that H' converges, hence remains bounded 
uniformly in S. 

Second, recall that G<jQ^{-,y) is the Green function in 9Qs with singularity at y. 
Since Rs equals on the boundary. Proposition |4.9| implies 



thus giving 



Therefore, 



Rs{x) = ^ ARsiy)G9Q,ix,y), (5.7) 



VRs{x) = ^Rs{y)yxGgQ,{x,y) 

Vi9Qs 



\vRs{x)\ = E I E ARs{y)VxGgQ,{x,y)\ 

< Y ^R^y) E \yxG,Q,{x,y)\ 

ye9Qs xeQs 

< Y ^R^y) E G9QA^,y) 

yi9Qs xeQs 

= Cs^ E E ARsiy)GsQAx,y) 

The second line uses the fact that ARs > 0, the third Proposition |4.10| the fifth Propo- 
sition 



4.9 



again, and the last inequality the facts that Qs contains of order 1/(5^ sites 
and that Rs is bounded uniformly in S (since Hs and Ss are). 

Thus, <5 ZajeQs |Vi?*| is uniformly bounded. Since the same result holds for 
{Fs/\/26)s>Q is precompact on Q (and more generally on any compact subset of fl) and 
the proof is completed. □ 



5.2 Convergence of the spin fermionic observable 

We now turn to the proof of convergence for the spin fermionic observable. Fix a simply 
connected domain {n,a,b) with two points on the boundary. For S > 0, always consider 
the spin fermionic observable on the discrete spin Dobrushin domain ($7^,05,65). Since 
the domain is fixed, we set Fs = Ffi<> as,bs- We follow the same three steps as before, 
beginning with the s-holomorphicity. The other two steps are only sketched, since they 
are more technical than in the FK-Ising case, see |CS09| . 

Proposition 5.7. For S > 0, Fs is s-holomorphic on fi^. 
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1(a) 1(6) 2 3(a) 3(6) 4 




Figure 13: The different possible cases in the proof of Proposition |5.7| w is depicted on 
the top, and uj' on the bottom. 

Proof Let x^y two adjacent medial vertices connected by the edge e = [xy]. Let v be 
the vertex of Q.s bordering the (medial) edge e. As before, set (resp. y^^) for the 
contribution of w to Fs{x) (resp. Fs{y)). We wish to prove that 

UJ UJ 

Note that the curve 7(0;) finishes at Xi^ or at so that uj cannot contribute to Fs{x) 
and Fs{y) at the same time. Thus, it is sufficient to partition the set of configurations 
into pairs of configurations (0;,^;'), one contributing to y, the other one to x, such that 

Pe(e){Xuj) = Pl{e){yuj')- 

Without loss of generality, assume that e is pointing southeast, thus £{e) = M (other 
cases can be done similarly). First note that 

z 

where 7(0^) is the interface in the configuration uj, 7' is any curve from as to bg (re- 
call that WY{as,bs) does not depend on 7'), and Z is a normalizing real number not 
depending on the configuration. There are six types of pairs that one can create, see 
Fig. [13] depicting the four main cases. Case 1 corresponds to the case where the interface 
reaches x or y and then extends by one step to reach the other vertex. In Case 2, 7 
reaches v before x and y, and makes an additional step to x or y. In Case 3, 7 reaches x 
or y and sees a loop preventing it from being extended to the other vertex (in contrast 
to Case 1). In Case 4, 7 reaches x or y, then goes away from v and comes back to the 
other vertex. Recall that the curve must always go to the left: in cases 1(a), 1(b), and 
2 there can be a loop or even the past of 7 passing through v. However, this does not 
change the computation. 

We obtain the following table for x^ and y^j' (we always express j/^/ in terms of 
x^). Moreover, one can compute the argument modulo tt of contributions Xi^ since the 
orientation of e is known. When upon projecting on M, the result follows. 
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configuration 


Case 1(a) 


Case 1(b) 


Case 2 


Case 3(a) 


Case 3(b) 


Case 4 


















(V2-l)e*"/V 


!jr/4 


e-'"/*a:„ 






e-5*"/*a;„ 


arg. Xu, mod n 


5n/S 


7r/8 


7r/8 


57r/8 


Sn/S 


57r/8 



□ 



Proof of Theorem 2.11 (Sketch). The proof is roughly sketched. We refer to |CS09| 
for a complete proof. 

Since Fs is s-harmonic, one can define the observable Hs as in Theorem |4.19| with 
the requirement that it is equal to on the white face adjacent to b. Then, Hg is constant 
equal to on the boundary as in the FK-Ising case. Note that Hs should not converge 
to 0, even if boundary conditions are away from a. Firstly, Hg is superharmonic and 
not harmonic, even though it is expected to be almost harmonic (away from a, H' and 
Hg are close), this will not be true near a. Actually, Hs should not remain bounded 
around a. 

The main difference compared to the previous section is indeed the unboundedness 
of Hs near as which prevents us from the immediate use of Proposition |4 . 3 [ It is actually 
possible to prove that away from a, Hs remains bounded, see |CS09| . This uses more 
sophisticated tools, among which are the boundary modification trick (see |DCHN10] 
for a quick description in the FK-Ising case, and [CS09 for the Ising original case). As 
before, boundedness implies precompactness (and thus boundedness) of {Fs)s>o away 
from a via Proposition |4.3[ Since Hs can be expressed in terms of Fs, it is easy to 
deduce that Hs is also precompact. 

Now consider a convergent subsequence {fs^iHs^) converging to {f,H). One can 
check that H is equal to on dil \ {a}. Moreover, the fact that Hg equals on the 
boundary and is superharmonic implies that Hg is greater than or equal to everywhere, 
implying H > m fl. This property of harmonic functions in a domain almost determines 
them. There is only a one-parameter family of positive harmonic functions equal to 
on the boundary. These functions are exactly the imaginary parts of conformal maps 
from n to the upper half-plane EI mapping a to oo. We can further assume that b is 
mapped to 0, since we are interested only in the imaginary part of these functions. 

Fix one conformal map ip from J7 to H, mapping a to oo and 6 to 0. There exists 
A > such that H = Xlmip. As in the case of the FK-Ising model, one can prove that 
Im (/^ f) = H, implying that f = \il>' . Since f{h) = 1 (it is obvious from the definition 

that Fs{bs) = 1), A equals . In conclusion, f{z) = \Jip'{z)/ip'{b) for every z eO.. □ 

Note that some regularity hypotheses on the boundary near b are needed to ensure 
that the sequence {fs^,Hs^) also converges near b. This is the reason for assuming that 
the boundary near b is smooth. We also mention that there is no normalization here. 
The normalization from the point of view of b was already present in the definition of 
the observable. 



6 Convergence to chordal SLE(3) and chordal SLE(16/3) 

The strategy to prove that a family of parametrized curves converges to SLE(k) follows 
three steps: 

• First, prove that the family of curves is tight. 
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• Then, show that any sub-sequential limit is a time-changed Loewner chain with a 
continuous driving process (see BefFara's course for details on Loewner chains and 
driving processes). 

• Finally, show that the only possible driving processes for the sub-sequential limits 
is \/~fiBt where Bt is a standard Brownian motion. 

The conceptual step is the third one. In order to identify the Brownian motion as 
being the only possible driving process for the curve, we find computable martingales 
expressed in terms of the limiting curve. These martingales will be the limits of fermionic 
observables. The fact that these (explicit) functions are martingales allows us to deduce 
martingale properties of the driving process. More precisely, we aim to use Levy's 
theorem: a continuous real-valued process X such that Xt and X"^ - at are martingales 
is necessarily -^aBt. 

6.1 Tightness of interfaces for the FK-Ising model 

In this section, we prove the following theorem: 

Theorem 6.1. Fix a domain (i},a,b). The family {'js)s>o of random, interfaces for 
the critical FK-Ising model in (f2,a,6) is tight for the topology associated to the curve 
distance. 

The question of tightness for curves in the plane has been studied in the groundbreak- 
ing paper |AB99| . In that paper, it is proved that a sufficient condition for tightness is 
the absence, at every scale, of annuli crossed back and forth an unbounded number of 
times. 

More precisely, for a: € and r < R, let Sr,R{x) = {x + [-R,R]^) \ {x + [-r, r]^) and 
define Ak{x;r, R) to be the event that there exist k crossings of the curve 75 between 
outer and inner boundaries of Sr,R{x). 

Theorem 6.2 (Aizenman-Burchard |AB99| ). Let ft be a simply connected domain and 
let a and b be two marked points on its boundary. Denote by ¥s the law of a random 
curve "fs on ilg from ag to bg. If there exist k eN, < 00 and A^, > 2 such that for all 
5 <r < R and xeVL, 

Vg{Ak{x-r,R))<Ck('^ \ 

then the family of curves (75) is tight. 

We now show how to exploit this theorem in order to prove Theorem |6.1| The main 
tool is Theorem 13. 161 

Lemma 6.3 (Circuits in annuli). Let £{x,n,N) be the probability that there exists an 
open path connecting the boundaries of Sn,N{x) . There exists a constant c< 1 such that 
for all n> 0, 

'^L,S„,2„(2;)(^(2^;"'27i)) < c. 

Note that the boundary conditions on the boundary of the annulus are wired. Via 
comparison between boundary conditions, this implies that the probability of an open 
path from the inner to the outer boundary is bounded uniformly on the configuration 
outside of the annulus. This uniform bound allows us to decouple what is happening 
inside the annulus with what is happening outside of it. 
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Figure 14: Left: The event v4g(a;, r, i?). In the case of exploration paths, it impHes the 
existence of alternating open and closed paths. Right: Rectangles Rt, Rr, Rb and Rj^ 
crossed by closed paths in the longer direction. The combination of these closed paths 
prevents the existence of a crossing from the inner to the outer boundary of the annulus. 



Proof Assume x = 0. The result follows from Theorem 3.16 (proved in Section 7.2) 



applied in the four rectangles Rb = [-2n,2n] x [-n, -2n], Rl = [-2n,-n] x [-2n,2n], 



Rt = [-2n, 2n] x [n, 2n] and Rr = [n, 2n] x [-2n, 2n], see Fig. 14 Indeed, if there exists 
a closed path crossing each of these rectangles in the longer direction, one can construct 
from them a closed circuit in Sn.2n- Now, consider any of these rectangles, Rb for 
instance. Its aspect ratio is 4, so that Theorem |3 . 1 6 1 implies that there is a closed path 
crossing in the longer direction with probability at least ci > (the wired boundary 



conditions are the dual of the free boundary conditions). The FKG inequality (3.2) 
implies that the probability of a circuit is larger than cf > 0. Therefore, the probability 
of a crossing is at most c=\- c\<\. □ 



We are now in a position to prove Theorem 6.1 



Proof of Theorem |6.1| Fix xeil, 5 < r < R and recall that we are on a lattice of 
mesh size 6. Let k to be fixed later. We first prove that 

<^^/;L(-^2fc(a;;r,2r))<c^- (6.1) 

for some constant c< 1 uniform in x,k,r,d and the configuration outside of Sr,2r{x). 

If A2k{^\ i"i 2r) holds, then there are (at least) k open paths, alternating with k dual 
paths, connecting the inner boundary of the annulus to its outer boundary. Since the 
paths are alternating, one can deduce that there are k open crossings, each one being 
surrounded by closed crossings. Hence, using successive conditionings and the compar- 
ison between boundary conditions, the probability for each crossing is smaller than the 
probability that there is a crossing in the annulus with wired boundary conditions (since 
these boundary conditions maximize the probability of £{x;r,2r)). We obtain 

'^a;L(-^2fc(a;;r,2r))<[0i,^,^^^^(,)(£:(x^ 

Using Lemma 6.3 4>\ ^ s 2 {x)(^(^'^ '^'')) ^ c < 1 and and ( |6.1[ ) follows. 



One can further fix k large enough so that c*^ < |. Now, one can decompose the 
annulus Sr,R{x) into roughly hi2{R/r) annuli of the form Sr,2r{x), so that for the 
previous k, 



^(^2fe(x;r,i?))<(-^) . (6.2) 



45 



Hence, Theorem 6.2 implies that the family (75) is tight. 



6.2 sub-sequential limits of FK-Ising interfaces are Loewner 
chains 

This subsection requires basic knowledge of Loewner chains and we refer to BefFara's 
course in this volume for an overview on the subject. In the previous subsection, 
traces of interfaces in Dobrushin domains were shown to be tight. The natural dis- 
crete parametrization does not lead to a suitable continuous parametrization. We would 
prefer our sub-sequential limits to be parametrized as Loewner chains. In other words, 
we would like to parametrize the curve by its so-called /i-capacity. In this case, we say 
that the curve is a time-changed Loewner chain. 

Theorem 6.4. Any sub-sequential limit of the family {'js)5>o of FK-Ising interfaces is 
a time-changed Loewner chain. 

Not every continuous curve is a time-changed Loewner chain. In the case of FK 
interfaces, the limiting curve is fractal-like and has many double points, so that the fol- 
lowing theorem is not a trivial statement. A general characterization for a parametrized 
non-selfcrossing curve in (Q,a,b) to be a time-changed Loewner chain is the following: 

• its /i-capacity must be continuous, 

• its /i-capacity must be strictly increasing. 

• the curve grows locally seen from infinity in the following sense: for any t > and 
for any e > 0, there exists 5 > such that for any s <t, the diameter of gs{ils'^^s+5) 
is smaller than e, where fls is the connected component of \ 7[0, s] containing b 
and gs is the conformal map from il^ to EI with hydrodynamical renormalization 
(see BefFara's course). 

The first condition is automatically satisfied by continuous curves. The third one usually 
follows from the two others when the curve is continuous, so that the crucial condition 
to check is the second one. This condition can be understood as being the fact that the 
tip of the curve is visible from b at every time. In other words, the family of hulls created 
by the curve (i.e. the complement of the connected component of f2 \ 74 containing b) 
is strictly increasing. This is the case if the curve does not enter long fjords created by 
its past at every scale, see Fig. [15] 

In the case of FK interfaces, this corresponds to so-called six arm event, and it boils 
down to proving that Ag > 2. A general belief in statistical physics is that many ex- 
ponents, called universal exponents, do not depend on the model. For instance, the 
so-called 5-arm exponent should equal 2. This would imply that Ag > A5 = 2. Unfor- 
tunately, computing the 5-arm exponent for the FK-Ising model is not an easy task. 
Therefore, we need to invoke a stronger structural theorem to prove that sub-sequential 
limits are Loewner chains. Recently, Kemppainen and the second author proved the 
required theorem, and we describe it now. 

For a family of parametrized curves {"fs)s>Oi define Condition (*) by: 

Condition {*): There exist C > 1 and A > such that for any < 5 < r < R/C , for 
any stopping time t and for any annulus Sr^R{x) not containing 7,-, the probability that 
75 crosses the annulus Sr^R{x) (from the outside to the inside) after tim e t while it is 



not forced to enter Sr.Ri^) again is smaller than C{r/R) , see Fig. 15 
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Figure 15: Left: An example of a fjord. Seen from b, the /i-capacity (roughly speakmg, 
the size) of the hull does not grow much while the curve is in the fjord. The event 
involves six alternating open and closed crossings of the annulus. Right: Conditionally 
on the beginning of the curve, the crossing of the annulus is unforced on the left, while 
it is forced on the right (it must go ultimately to b). 



Roughly speaking, the previous condition is a uniform bound on unforced crossings. 
Note that it is necessary to assume the fact that the crossing is unforced. 



Theorem 6.5 f |KS10| l. // a family of curves (75) satisfies Condition {*), then it is 
tight for the topology associated to the curve distance. Moreover, any sub-sequential limit 
(75^) is to a time- changed Loewner chain. 

Tightness is almost obvious, since Condition (*) implies the hypothesis in Aizenman- 
Burchard's theorem. The hard part is the proof that Condition (*) guarantees that the 
/i-capacity of sub-sequential limits is strictly increasing and that they create Loewner 
chains. The reader is referred to [KSIO] for a proof of this statement. We are now in a 
position to prove Theorem |6.4[ 



Proof of Theorem 6.4 Lemma 6.3 allows us to prove Condition (*) without difficulty. 
□ 



6.3 Convergence of FK-Ising interfaces to SLE(16/3) 

The FK fermionic observable is now proved to be a martingale for the discrete curves 
and to identify the driving process of any sub-sequential limit of FK-Ising interfaces. 

Lemma 6.6. Let 5 > 0. The FK fermionic observable A/^(z) = ^o^x-y[o,n],7„,b5 (-2^) ^ 
martingale with respect to (J-n), where Tn is the a-algebra generated by the FK interface 
7[0,n]. 

Proof For a Dobrushin domain (il^, 05,65), the slit domain created by "removing" 
the first n steps of the exploration path is again a Dobrushin domain. Conditionally on 
7[0,n], the law of the FK-Ising model in this new domain is exactly 0o"'''^rn t This 

observation implies that M^{z) is the random variable l^g^^e^*^^^^^'''' conditionally on 
J-n, therefore it is automatically a martingale. □ 



Proposition 6.7. Any sub-sequential limit of {'js)s>o which is a Loewner chain is the 
(chordal) Schramm- Loewner Evolution with parameter k = 16/3. 
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Proof Consider a sub-sequential limit 7 in the domain {n,a,b) which is a Loewner 
chain. Let be a map from (fl, a, b) to (H, 0, 00). Our goal is to prove that 7 = 0(7) is 
a chordal SLE(16/3) in the upper half-plane. 

Since 7 is assumed to be a Loewner chain, 7 is a growing hull from to 00 
parametrized by its /i-capacity. Let Wt be its continuous driving process. Also, de- 
fine gt to be the conformal map from IHI\7[0,<] to H such that gt{z) = z + 2t/z + 0{l/z^) 
when z goes to 00. 

Fix z' € fl. For S > 0, recall that M^(z') is a martingale for 75. Since the martingale 
is bounded, M^^{z') is a martingale with respect to J^n, where Tt is the first time 
at which (jjils) has an /i-capacity larger than t. Since the convergence is uniform, 
Mt{z') := \ims->o M^^{z') is a martingale with respect to Qt, where Gt is the cr-algebra 
generated by the curve 7 up to the first time its /i-capacity exceeds t. By definition, this 
time is t, and Gt is the a-algebra generated by 7[0,t]. 

Recall that Mt{z') is related to (l>{z') via the conformal map from H \ 7[0,i] to 
Mx(0, 1), normalized to send jt to -00 and 00 to 00. This last map is exactly ^ ln.{gt-Wt)- 
Setting z = (j){z'), we obtain that 



V^M[- := V^Mtiz') = ^[\nigtiz)-Wt)] 



9'ti^) 



is a martingale. Recall that, when z goes to infinity, 

2t 



9t{z) 



z + 



^^ 9t{z)-wt 

^.0(1) and 



(6.3) 



(6.4) 



For s < t. 



^■¥.[M^\Gs'\ = E 



Taking s = t yields 



l-2t/z2+0(l/23) 

\| z-Wt + 2t/z + 0{l/z^) 



E 



1 + Wt/z + 1 {3W^ - m) /z^ + O (l/z^) I Gs 



(1 + ^E[Wt\GsVz + Ie[3w^ - miGsVz^ + O (l/z^)) 

^ (1 + \Ws/z + ^{SW! - 16s)/z' + 0(l/z^)) . 



Since E[M(^|5s] = M|, terms in the previous asymptotic development can be matched 
together so that E[M^t|e;,] = Ws and E[W^ - ^t\Gs] = - ^s. Since Wt is continuous. 

Levy's theorem implies that Wt = \J^Bt where Bt is a standard Brownian motion. 

In conclusion, 7 is the image by (p~^ of the chordal Schramm-Loewner Evolution with 
parameter k = 16/3 in the upper half-plane. This is exactly the definition of the chordal 
Schramm-Loewner Evolution with parameter k = 16/3 in the domain (f2,a,6). □ 



Proof of Theorem |3.13| By Theorem |4.3| the family of curves is tight. Using The- 
any sub-sequential limit is a time-changed Loewner chain. 



orem 



6.4 



Consider such a 

sub-sequential limit and parametrize it by its /i-capacity. Proposition |6.7| then implies 
that it is the Schramm-Loewner Evolution with parameter n = 16/3. The possible limit 
being unique, the claim is proved. □ 
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6.4 Convergence to SLE(3) for spin Ising interfaces 

The proof of Theorem [2l0| is very similar to the proof of Theorem |3.13| except that 
we work with the spin Ising fermionic observable instead of the FK-Ising model one. 
The only point differing from the previous section is the proof that the spin fermionic 
observable is a martingale for the curve. We prove this fact now and leave the remainder 
of the proof as an exercise. Let 7 be the interface in the critical Ising model with 
Dobrushin boundary conditions. 

Lemma 6.8. Let 6>0, the spin fermionic observable M^(z) = i^o»N-y[o.n].7(n),b5 (-z) a 
martingale with respect to {Tn), where Tn is the a-algebra generated by the exploration 
process 7[0, n]. 



Proof It is sufficient to check that Fs{z) has the martingale property when 7 = 'y{uj) 
makes one step 71. In this case J^o is the trivial cr-algebra, so that we wish to prove 

^^l^fi [-^f^Mii7i],7i.fc5(^)] = Fn'iMsMi^)' (6-5) 

where /^^'''f^ is the critical Ising measure with Dobrushin boundary conditions in fl. 
Write Zqo^ i^^ i,^ (resp. Z^o^^g^^^^ ^^ i,^) for the partition function of the Ising model 
with Dobrushin boundary conditions on {n'^,as,bs) (resp. (17* \ [asx],x,bs)), i.e. 
Zn'>^[asx],x,bs = Et^(\/2- l)'"^'. Note that Zo«x[a,a;],x,6a is almost the denominator of 
Fni^[asx],xMi^s)- By definition, 

Zn^.as.bs M^l'^n (71 =X)= (\/2 - l)Zi'io^[asx],x,bs 



FQ.l^Ya^x\,xf,f,{.Z&) 



: e'2 



In the second equality, we used the fact that f n°\[a5x] (2^1 ^5) is in bijection with config- 
urations of £^Qo (05,25) such that 71 = a; (there is still a difference of weig ht of \/2 - 1 
between two associated configurations). This gives 



t^0^ni^i = ^)Fn^,-[asx],x.bAzs} = -q^^r^^;^^ ■ 



The same holds for all possible first steps. Summing over all possibilities, we obtain the 
expectation on one side of the equality and Fqo ,^^ i,^{zs) on the other side, thus proving 

" ' a 



Exercise 6.9. Prove that spin Ising interfaces converge to SLE(3). For tightness and 
the fact that sub-sequential limits are Loewner chains, it is sufficient to check Condition 
(*). To do so, try to use Theorem \ 3.1()\ and the Edwards-Sokal coupling to prove an 
intermediate result similar to Lemma \6. M 



7 Other results on the Ising and FK-Ising models 
7.1 Massive harmonicity away from criticality 

In this subsection, we consider the fermionic observable F for the FK-Ising model away 
from criticality. The Ising model is still solvable and the observable becomes massive 
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harmonic (i.e. A/ = A^/). We refer to |BDC11| for details on this paragraph. We start 
with a lemma which extends Lemma 5.2 to p psd = \/2/(l + \/2). 



Lemma 7.1. Let p e (0, 1) . Consider a vertex v e ft^ \ dil" , 

F{A)-F{C) = ic"' [F{B) - F{D)] (7.1) 

where A is an adjacent (to v) medial edge pointing towards v and B, C and D are 
indexed in such a way that A, B, C and D are found in counterclockwise order. The 
parameter a is defined by 

^,„_ e-'-/4(l-p)x/2+p 

The proof of this statement follows along the same lines as the proof of Lemma |5.2[ 

Proposition 7.2. For p < \/2/(l + \/2), there exists ^ = ^(p) > such that for every n, 

(j)p{0^m) < e"«", (7.2) 

where the mesh size of the lattice h is 1. 

In this proof, the lattices are rotated by an angle tt/4. We will be able to estimate 
the connectivity probabilities using the FK fermionic observable. Indeed, the observable 
on the free boundary is related to the probability that sites are connected to the wired 
arc. More precisely: 

Lemma 7.3. Fix (G,a, 5) a Dobrushin domain and p e (0, 1). Let u e G be a site on 
the free arc, and e be a side of the black diamond associated to u which borders a white 
diamond of the free arc. Then, 

\F{e)\ = (j)'^'Q{u ^ wired a,ic). (7.3) 

Proof Let u be a site of the free arc and recall that the exploration path is the interface 
between the open cluster connected to the wired arc and the dual open cluster connected 
to the free arc. Since u belongs to the free arc. u is connected to the wired arc if and 
only if e is on the exploration path, so that 

The edge e being on the boundary, the exploration path cannot wind around it, so that 
the winding (denoted Wi) of the curve is deterministic (and easy to write in terms of 
that of the boundary itself). We deduce from this remark that 

\F{e)\ = |0;;^(e4Wu,,,)| = {e^"^^ <t>;%{e e j)\ 

= <^p,g(s ^ t) = <Pp!g{u ^ wired arc). 

□ 

We are now in a position to prove Proposition |7.2| We first prove exponential decay 
in a strip with Dobrushin boundary conditions, using the observable. Then, we use 
classical arguments of FK percolation to deduce exponential decay in the bulk. We 
present the proof quickly (see [BDCll^ for a complete proof). 
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Proof Let p < psd, and consider the FK-Ising model of parameter p in the strip of 
height £, with free boundary conditions on the top and wired boundary conditions on 
the bottom (the measure is denoted by (t'^g^)- It is easy to check that one can define 
the FK fermionic observable F in this case, by using the unique interface from -oo to 
oo. This observable is the limit of finite volume observables, therefore it also satisfies 
Lemma 17.11 

Let efe be the medial edge with center ik + see Fig. 
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A simple computation 

using Lemmata 7.1 and 5.1 plus symmetries of the strip (via translation and horizontal 

[1 + cos(7r/4 - a)] cos(7r/4 - a) 



reflection) implies 



[1 + cos(7r/4 + a)] cos(7r/4 + a) 
Using the previous equality inductively, we find for every i > 

\F{ee)\ = e-^'\F{eo)\<e-i' 



(7.4) 



with 



■In 



[1 + cos(7r/4 - a)] cos(7r/4 - a) 



[1 + cos(7r/4 + a)] cos(7r/4 + a) ' 
Since ei is adjacent to the free arc. Lemma |7. 3| implies 

; oo , — I 



(7.5) 



•[i£-Z] = |F(e,)|<e-«^ 



a = (01,02) 




Figure 16: Left: Edges and Ck+i- Right: A dual circuit surrounding an open 
path in the box [-02,02]^. Conditioning on to the most exterior such circuit gives no 
information on the state of the edges inside it. 

Now, let iV € N and recall that (f>'^ j^i '■= 4'p2 [-n n]^ converges to the infinite-volume 
measure with free boundary conditions 0^ when N goes to infinity. 

Consider a configuration in the box [-N, N]'^, and let A^ax be the site of the cluster 
of the origin which maximizes the £°°-norm max{|a;i|, \x2\} (it could be equal to N). If 
there is more than one such site, we consider the greatest one in lexicographical order. 
Assume that Amax equals a = ai + \a2 with 02 > |ai| (the other cases can be treated the 
same way by symmetry, using the rotational invariance of the lattice). 

By definition, if Amax equals a, a is connected to in [-02,02]^. In addition to this, 
because of our choice of the free boundary conditions, there exists a dual circuit starting 
from a + i/2 in the dual of [-02, a2]^ (which is the same as L* n [-02 - 1/2, 02 + 1/2]^) 
and surrounding both a and 0. Let F be the outermost such dual circuit: we get 

<jv(^max = a) = E <iv(« ^ = 7)0",w(r = 7), (7.6) 

7 
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where the sum is over contours 7 in the dual of [-02,02]^ that surround both a and 0. 

The event {F = 7} is measurable in terms of edges outside or on 7. In addition, 
conditioning on this event implies that the edges of 7 are dual-open. Therefore, from 
the domain Markov property, the conditional distribution of the configuration inside 7 is 
a FK percolation model with free boundary conditions. Comparison between boundary 
conditions implies that the probability of {a 0} conditionally on {F = 7} is smaller 
than the probability of {a 0} in the strip with free boundary conditions on the 
top and wired boundary conditions on the bottom. Hence, for any such 7, we get 

<l>%{a - 0|F = 7) < CC(" ^ 0) = - Z) < e-^"^^ 

(observe that for the second measure, Z is wired, so that {a ^ 0} and {a ^ Z} have the 
same probability). Plugging this into (7.6), we obtain 

0°,jv(^max = a) < ^e-«'"^'^{'^-'^^} ^%{r = 7) < e-^"^^ = e'^'^-^^-^'-^l 
7 

Fix n < N. We deduce from the previous inequality that there exists a constant 
< c < 00 such that 

0°,jv(O - Z2 X [-n, nf) < <l>pMAn.. = a)< cne'^^ . 

ae[-N,Ny\[-n,nY 

Since the estimate is uniform in N, we deduce that 

0°(O ^ in) < (/)°(0 ^ Z^ X [-n,nf) < cne-^'^ . (7.7) 

□ 



Theorem 7.4. The critical parameter for the FK-Ising model is \/2/(l + \/2). The 
critical inverse-temperature for the Ising model is | ln(l + \/2). 



Proof The inequality Pc > \/2/(l + \/2) follows from Proposition 7.2 since there is 
no infinite cluster for (p^ 2 when p < psd (the probability that and in are connected 
converges to 0). In order to prove that pc < \/2/(l + \/2), we harness the following 
standard reasoning. 

Let An be the event that the point n e N is in an open circuit which surrounds the 
origin. Notice that this event is included in the event that the point n e N is in a cluster 



of radius larger than n. For p < \/2/(l + %/2), a modification of (7.7 1 implies that the 
probability of An decays exponentially fast. The Borel-Cantelli lemma shows that there 
is almost surely a finite number of n such that An occurs. In other words, there is a.s. 
only a finite number of open circuits surrounding the origin, which enforces the existence 
of an infinite dual cluster whenever p < \/2/(l + \/2). Using duality, the primal model is 
supercritical whenever p > \/2/ ( 1 + \/2) , which implies p^ < \/2/ ( 1 + \/2) . □ 

In fact, the FK fermionic observable Fg in a Dobrushin domain {flg,as,bs) is massive 
harmonic when p + psd- More precisely. 

Proposition 7.5. Letp + psd, 

AsFs{v) = {cos2a-l)Fs{v) (7.8) 

for every v e fl"^ \dfl'^, where Ag is the average on sites at distance \/25 minus the value 
at the point. 
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When 5 goes to 0, one can perform two scaling limits. If p = psd{^-^S) goes to psd as 
S goes to 0, j2 (^5 + [1 - cos 2a] J) converges to A + X^I. Then Fs (properly normalized) 
should converge to a function / satisfying A/ + X^f = inside the domain. Except for 
A = 0, the limit will not be holomorphic, but massive harmonic. Discrete curves should 
converge to a limit which is not absolutely continuous with respect to SLE(16/3). The 
study of this regime, connected to massive SLEs, is a very interesting subject. 

If we fix p < psci, one can interpret massive harmonicity in terms of killed random 
walks. Roughly speaking, Fs{v) is the probability that a killed random walk starting at 
V visits the wired arc dba- Large deviation estimates on random walks allow to compute 
the asymptotic of Fg inside the domain. In |BDC11| . a surprising link (first noticed by 
Messikh jMesOH]) between correlation lengths of the Ising model and large deviations 
estimates of random walks is presented. We state the result in the following theorem: 

Theorem 7.6. Fix f3 < Pc (and a associated to it) and set 

m{f5) := cos(2a). 



For any a; e L, 



- lim —h\iip\aQ(j{nxy\ = - lim — In G'm(^)(0, nx). (7.9) 



Above, Gm{0,x) := E^[to'^] for any x € L and m < 1, where r is the hitting time of the 
origin and is the law of a simple random walk starting at x. 



The massive Green function 6*^(0,2;) on the right of (7.9 1 has been widely studied 



In particular, we can compute the rate of decay in any direction and deduce Theorem|2.3 



and Theorem [0](see e.g. |Mes06| ). 



Exercise 7.7. Prove Lemma \ 7.1\ and the fact that F is massive harmonic inside the 
domain. 

7.2 Russo-Seymour- Welsh Theorem for FK-Ising 



In this section, we sketch the proof of Theorem 3.16 see [DCHNIO] for details. We 
would like to emphasize that this result does not make use of scaling limits. Therefore, 
it is mostly independent of Sections [5] and |6] 

We start by presenting a link between discrete harmonic measures and the probability 
for a point on the free arc dab of a FK Dobrushin domain to be connected to the wired 
arc dba- 

Let us first define a notion of discrete harmonic measure in a FK Dobrushin domain 
which is slightly different from the usual one. First extend Vlg u fij by adding two 
extra layers of vertices: one layer of white faces adjacent to 9*^, and one layer of black 
faces adjacent to dba- We denote the extended domains by (is and Vt*g. 

Define (X*)j>q to be the continuous-time random walk on the black faces that jumps 
with rate 1 on neighbors, except for the faces on the extra layer adjacent to dab onto 
which it jumps with rate p •■= 2/(\/2+l). For B e fi^, we denote by H'{B) the probability 
that the random walk X' starting at B hits dil.s on the wired arc dba (in other words, if 
the random walk hits dba before hitting the extra layer adjacent to d*i^). This quantity 
is called the (modified) harmonic measure of dba seen from B. Similarly, one can define 
a modified random walk and the associated harmonic measure of 9*^, seen from w. 
We denote it by iJ°(w). 

Proposition 7.8. Consider a FK Dobrushin domain (fi^, 05,65), for any site B on the 
free arc dab, 

\/h"{W) < [B - dba] < \/h'{B), (7.10) 
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where W is any dual neighbor of B not on d*f^ . 

This proposition raises a connection between harmonic measure and connectivity 
properties of the FK-Ising model. To study connectivity probabihties for the FK-Ising 
model, it suffices to estimate events for simple random walks (slightly repelled on the 
boundary). The proof makes use of a variant of the "boundary modification trick". This 



trick was introduced in fCS09] to prove Theorem 2.11 It can be summarized as follows: 
one can extend the function iJ by or 1 on the two extra layers, then H' (resp. H°) 
is subharmonic (resp. superharmonic) for the Laplacian associated to the random walk 
X' (resp. X°). Interestingly, H' is not subharmonic for the usual Laplacian. This trick 
allows us to fix the boundary conditions (0 or 1), at the cost of a slightly modified notion 
of harmonicity. 



We can now give the idea of the proof of Theorem 3.16| (we refer to [DCHNIO] for 



details) . The proof is a second moment estimate on the number of pairs of connected 



sites on opposite edges of the rectangle. We mention that another road to Theorem 3.16 
has been proposed in |KS10| . 



Proof of Theorem 3.16 (Sketch) Let i?„ = [0,4n] x [0,7i] be a rectangle and let 
N be the number of pairs {x,y), x 6 {0} x [0,n] and y 6 {4n} x [0,n] such that x 
is connected to y by an open path. The expectation of N is easy to obtain using the 
previous proposition. Indeed, it is the sum over all pairs x, y of the probability of {x ^ y} 
when the boundary conditions are free. Free boundary conditions can be thought of as 
a degenerate case of a Dobrushin domain, where a = b = y. In other words, we want to 
estimate the probability that x is connected to the wired arc d^a = {y}- Except when x 
and y are close to the corners, the harmonic measure of y seen from x is of order 
so that the probability of {x <-> y} is of order 1/n. Therefore, there exists a universal 
constant c > such that 'Pp ^ n [^] ^ cn. 

The second moment estimate is harder to obtain, as u sual. Nevertheless, it can be 



proved, using successive conditioning and Proposition 7.8 that (t^p ^ r [-^1 - ^"-"^ 
some universal C > 0; see [DCHNlO] for a complete proof. Using the Cauchy-Schwarz 
inequality, we find 



which implies 



uniformly in n. 



4^^_p^j3 open crossing] = 4>R„.p^^[N > Q] > c^/C (7.12) 



We have already seen that Theorem 3.16 is central for proving tightness of interfaces. 



We would also like to mention an elementary consequence of Theorem |3.16| 
Proposition 7.9. There exist constants < c, C, (5, A < oo such that for any sites y e L, 

where fip^ is the unique infinite-volume measure at criticality. 



Proof Using the Edwards-Sokal coupling, (7.131 can be rephrased as 



4f ' - 
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where 4'psd,2 is the unique FK-Ising infinite- volume measure at criticahty. In order to 
get the upper bound, it suffices to prove that (/)p^^,2(0 dAk) decays polynomially 
fast, where Ak is the box of size k = \x - y\ centered at x. We consider the annuli 
An = S'2"-i.2" (a;) for n < ln2 k, and £{An) the event that there is an open path crossing 



An from the inner to the outer boundary. We know from Corollary 6.3 (which is a direct 



application of Theorem 3.161 that there exists a constant c < 1 such that 

</'A„,P..,2(^(^n))^C 

for all n > 1. By successive conditionings, we then obtain 

ln2 k 
n=l 

and the desired result follows. The lower bound can be done following the same kind of 
arguments (we leave it as an exercise). □ 

Therefore, the behavior at criticality (power law decay of correlations) is very differ- 
ent from the subcritical phase (exponential decay of correlations). Actually, the previous 
result is far from optimal. One can compute correlations between spins of a domain very 
explicitly. In particular, /i^^ [fja-cry] behaves like jx-yl"", where a = 1/4. We mention 
that a is one example of critical exponent. Even though we did not discuss how compute 
critical exponents, we mention that the technology developed in these notes has for its 
main purpose their computation. 

To conclude this section, we mention that Theorem |3.16| leads to ratio mixing prop- 



erties (see Exercise 7.101 of the Ising model. Recently, Lubetzky and Sly |LS10| used 
these spatial mixing properties in order to prove an important conjecture on the mixing 
time of the Glauber dynamics of the Ising model at criticality. 

Exercise 7.10 (Spatial mixing). Prove that there exist c, A > such that for any r < R, 

|0p^„2(AnB)-0p^,,2(A)0p^,^2(S)|<c(^)'^0p.,.2(A)0p^„2(S) (7.14) 
for any event A (resp. B) depending only on the edges in the box [-r, r]^ (resp. outside 

7.3 Discrete singularities and energy density of the Ising model 

In this subsection, we would like to emphasize the fact that slight modifications of 
the spin fermionic observable can be used directly to compute interesting quantities of 
the model. Now, we briefiy present the example of the energy density between two 



neighboring sites x and y (Theorem 2.121 



So far, we considered observables depending on a point a on the boundary of a 
domain, but we could allow more ffexibility and move a inside the domain: for e fi^, 
we define the fermionic observable F^^^^{zs) for t as by 

F^^izs) = A ' (7.15) 

where A is a well-chosen explicit complex number. Note that the denominator of the 
observable is simply the partition function for free boundary conditions g. Actually, 
using the high-temperature expansion of the Ising model and local rearrangements, the 
observable can be related to spin correlations |HS10| : 
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Lemma 7.11. Let [xy] be an horizontal edge of fig- Then 



^ Mi,n,K^J.] = Pi{ac)[F,l (C)] + Piiad)[F^, (d)] 

where a is the center of [xy], c = a + and d = a - 5^^. 

If A is chosen carefully, the function F^^ is s-holoniorphic on Q,s n {as}. Moreover, 
its complex argument is fixed on the boundary of the domain. Yet. the function is 
not s-holomorphic at as (meaning that there is no way of defining i^^ {as) so that the 
function is s-holomorphic at as). In other words, there is a discrete singularity at a, 
whose behavior is related to the spin-correlation. 

We briefly explain how one can address the problem of discrete singularities, and 
we refer to |HS10| for a complete study of this case. In the continuum, singularities 
are removed by subtracting Green functions. In the discrete context, we will do the 
same. We thus need to construct a discrete s-holomorphic Green function. Preholo- 
morphic Green function^ were already constructed in [KenOO' . These functions are not 
s-holomorphic but relevant linear combinations of them are, see [HSIO . We mention 
that the s-holomorphic Green functions are very explicit and their convergence when 
the mesh size goes to can be studied. 



Proof of Theorem 2.12 (Sketch) The function F^^ jS converges uniformly on any 
compact subset of VL \{a|. This fact is not helpful, since the interesting values of F^'' are 
located at neighbors of the singularity. It can be proved that, subtracting a well-chosen 
s-holomorphic Green function g"^^, one can erase the singularity at as. More precisely, 
one can show that [-F^* ~5o''^]/'^ converges uniformly on towards an explicit conformal 
map. The value of this map at a is ^^^(a). Now, /i^ uJl'^^'^vI '^^'^ expressed in 
terms of F^?* for neighboring vertices of ag. Moreover, values of g"^^ for neighbors of ag 
can be computed explicitly. Using the fact that 

Fr=9z+5-\[Fr-9zi 

and Lemma |7.11| the convergence result described above translates into the following 
asymptotics for the spin correlation of two neighbors 



8 Many questions and a few answers 
8.1 Universality of the Ising model 

Until now, we considered only the square lattice Ising model. Nevertheless, normaliza- 
tion group theory predicts that the scaling limit should be universal. In other words, the 
limit of critical Ising models on planar graphs should always be the same. In particular, 
the scaling limit of interfaces in spin Dobrushin domains should converge to SLE(3). 

Of course, one should be careful about the way the graph is drawn in the plane. 
For instance, the isotropic spin Ising model of Section |2] when considered on a stretched 
square lattice (every square is replaced by a rectangle), is not conformally invariant (it is 

®i.e. satisfying the Cauchy-Riemann equation except at a certain point. 
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not invariant under rotations) . Isoradial graphs form a large family of graphs possessing 
a natural embedding on which a critical Ising model is expected to be conformally 
invariant. More details are now provided about this fact. 

Definition 8.1. A rhombic embedding of a graph G is a planar quadrangulation satis- 
fying the following properties: 

• the vertices of the quadrangulation are the vertices of G and G* , 

• the edges connect vertices of G to vertices of G* corresponding to adjacent faces 
ofG, 

• all the edges of the quadrangulation have equal length, see Fig. | jQ[ 

A graph which admits a rhombic embedding is called isoradial. 

Isoradial graphs are fundamental for two reasons. First, discrete complex analy- 
sis on isoradial graphs was extensively studied (see e.g. [MerOll IKen02i ICS08| ) as 
explained in Section |4] Second, the Ising model on isoradial graphs satisfies very spe- 
cific integrability properties and a natural critical point can be defined as follows. Let 
Jxy = arctanh[tan (0/2)] where 6 is the half-angle at the corner x (or equivalently y) 
made by the rhombus associated to the edge [xy]. One can define the critical Ising 
model with Hamiltonian 

= - JxyCxCTy 

x~y 

This Ising model on isoradial graphs (with rhombic embedding) is critical and confor- 
mally invariant in the following sense: 

Theorem 8.2 (Chelkak, Smirnov |CS09| ). The interfaces of the critical Ising model on 
isoradial graphs converge, as the mesh size goes to 0, to the chordal Schramm- Loewner 
Evolution with k = 3. 

Note that the previous theorem is uniform on any rhombic graph discretizing a given 
domain (fi, a, 6), as soon as the edge-length of rhombi is small enough. This provides a 
first step towards universality for the Ising model. 

Question 8.3. Since not every topological quadrangulation admits a rhombic embed- 
ding IKS05I, can another embedding with a sufficiently nice version of discrete complex 
analysis always be found? 

Question 8.4. Is there a more general discrete setup where one can get similar esti- 
mates, in particular convergence of preholomorphic functions to the holomorphic ones 
in the scaling limit? 

In another direction, consider a biperiodic lattice C (one can think of the universal 
cover of a finite graph on the torus) , and define a Hamiltonian with periodic correlations 
{Jxy) by setting 11(a) = - Y.x-y JxyO'x'^y The Ising model with this Hamiltonian makes 
perfect sense and there exists a critical inverse temperature separating the disordered 
phase from the ordered phase. 

Question 8.5. Prove that there always exists an embedding of C such that the Ising 
model on C is conformally invariant. 
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cluster-weight q 



subcritical phase 



critical phase: first order 




edge-weight p 



Figure 17: The phase diagram of the FK percolation model on the square lattice. 



8.2 Full scaling limit of critical Ising model 

It has been proved in |KSrO| that the scaling limit of Ising interfaces in Dobrushin 
domains is SLE(3). The next question is to understand the full scaling limit of the 
interfaces. This question raises interesting technical problems. Consider the Ising model 
with free boundary conditions. Interfaces now form a family of loops. By consistency, 
each loop should look like a SLE(3). In jHKll| . Hongler and Kytola made one step 
towards the complete picture by studying interfaces with +f - /free boundary conditions. 

Sheffield and Werner |SW10a[ ISWlOb] introduced a one-parameter family of pro- 
cesses of non-intersecting loops which are conformally invariant - called the Conformal 
Loop Ensembles CLE(k) for k > 8/3. Not surprisingly, loops of CLE(k) are locally simi- 
lar to SLE(/c), and these processes are natural candidates for the scaling limits of planar 
models of statistical physics. In the case of the Ising model, the limits of interfaces all 
together should be a CLE(3). 

8.3 FK percolation for general cluster-weight g > 

The FK percolation with cluster-weight q e (0, oo) is conjectured to be critical for Pc{q) = 
y/q/{l + y/q) (see [BDClOj for the case q > I). Critical FK percolation is expected to 
exhibit a very rich phase transition, whose properties depend strongly on the value of 
q (see Fig. [Tt]). We use generalizations of the FK fermionic observable to predict the 
critical behavior for general q. 
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8.3.1 Case < q < 4. 

The critical FK percolation in Dobrushin domains can be associated to a loop model 
exactly like the FK-Ising model: each loop receives a weight ^Jq. In this context, one 
can define a natural generalization of the fermionic observable on medial edges, called 
a parafermionic observable, by the formula 

F{e) = E^,»,,,,,,p,,K-''^-(^'''^)le.^], (8.1) 



where a = a{q) is called the spin (a takes a special value described below). Lemma |5.2| 
has a natural generalization to any q e [0, oo): 

Proposition 8.6. For q < 4 and any FK Dobrushin domain, consider the observable F 
at criticality with spin a = 1- ^ arccos(^/q/2) . For any medial vertex inside the domain, 

F{N)-F{S) = i[F{E)- F{W)] (8.2) 

where N , E, S and W are the four medial edges adjacent to the vertex. 

These relations can be understood as Cauchy-Riemann equations around some ver- 
tices. Importantly, F is not determined by these relations for general q (the number of 
variables exceeds the number of equations). For q = 2, which corresponds to cr = 1/2, the 



complex argument modulo tt of the observable offers additional relations (Lemma 5.1) 



and it is then possible to obtain the preholomophicity (Proposition 5.3 1 



Parafermionic observables can be defined on medial vertices by the formula 

^ e~v 

where the summation is over medial edges with u as an endpoint. Even though they are 
only weakly-holomorphic, one still expects them to converge to a holomorphic function. 
The natural candidate for the limit is not hard to find: 

Conjecture 8.7. Let q < 4 and {Q,a,b) be a simply connected domain with two points 
on its boundary. For every z e fl, 

:Fs{z) -* (j)'{zY when 6 (8.3) 



m 

where a =1- - arccos(^yq/2), Fg is the observable (oi Pc{<i)) in discrete domains with 
spin a, and (p is any conformal map from to Rx (0, 1) sending a to -oo and b to oo. 

Being mainly interested in the convergence of interfaces, one could try to follow the 
same program as in Section [6] 

• Prove compactness of the interfaces. 

• Show that sub-sequential limits are Loewner chains (with unknown random driving 
process Wt). 

• Prove the convergence of discrete observables (more precisely martingales) of the 
model. 

• Extract from the limit of these observables enough information to evaluate the 
conditional expectation and quadratic variation of increments of Wt (in order to 
harness the Levy theorem). This would imply that Wt is the Brownian motion 
with a particular speed k and so curves converge to SLE(k). 
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The third step, corresponding to Conjecture |8.7[ should be the most difficult. Note 
that the ffist two steps are also open for q + 0,1,2. Even though the convergence of 
observables is still unproved, one can perform a computation similar to the proof of 



Proposition 6.7 in order to identify the possible limiting curves (this is the fourth step). 



The following conjecture is thus obtained: 

Conjecture 8.8. For g < 4, the law of critical FK interfaces converges to the Schramm- 
Loewner Evolution with parameter k = At:/ arccos(-^/g/2) . 

The conjecture was proved by Lawler, Schramm and Werner |LSW04a| for q = 0, 
when they showed that the perimeter curve of the uniform spanning tree converges to 
SLE(8). Note that the loop representation with Dobrushin boundary conditions still 
makes sense for q = (more precisely for the model obtained by letting q ^ and 
p/q ^ 0). In fact, configurations have no loops, just a curve running from a to b (which 
then necessarily passes through all the edges), with all configurations being equally 
probable. The q = 2 case corresponds to Theorem |3.13[ All other cases are wide open. 
The q = 1 case is particularly interesting, since it is actually bond percolation on the 
square lattice. 



8.3.2 Case g > 4. 

The picture is very different and no conformal invariance is expected to hold. The phase 
transition is conjectured to be of first order : there are multiple infinite- volume measures 
at criticality. In particular, the critical FK percolation with wired boundary conditions 
should possess an infinite cluster almost surely while the critical FK percolation with 
free boundary conditions should not (in this case, the connectivity probabilities should 
even decay exponentially fast). This result is known only for q > 25.72 (see |Gri06| and 
references therein). 

Note that the observable still makes sense in the g > 4 case, providing a is chosen 
so that 2sin(7r(j/2) = ^g. Interestingly, a becomes purely imaginary in this case. A 
natural question is to relate this change of behavior for a with the transition between 
conformally invariant critical behavior and first order critical behavior. 

8.4 0{n) models on the hexagonal lattice 

The Ising fermionic observable was introduced in [Smi06^ in the setting of general 0{n) 
models on the hexagonal lattice. This model, introduced in [D MNS81| on the hexagonal 
lattice, is a lattice gas of non-intersecting loops. More precisely, consider configurations 
of non-intersecting simple loops on a finite subgraph of the hexagonal lattice and intro- 
duce two parameters: a loop-weight n > (in fact n > -2) and an edge-weight x > 0, and 
ask the probability of a configuration to be proportional to 'o^p^^x* "^gcs 

Alternatively, an interface between two boundary points could be added: in this 
case configurations are composed of non-intersecting simple loops and one self-avoiding 
interface (avoiding all the loops) from a to b. 

The 0(0) model is the self-avoiding walk, since no loop is allowed (there is still a self- 
avoiding path from a to b). The 0(1) model is the high-temperature expansion of the 
Ising model on the hexagonal lattice. For integers n, the 0(n)-model is an approximation 
of the high-temperature expansion of spin 0(n)-models (models for which spins are n- 
dimensional unit vectors). 

The physicist Bernard Nienhuis |Nie82l lNie84| conjectured that 0(n)-models in the 
range n e (0, 2) (after certain modifications n € (-2, 2) would work) exhibit a Berezinsky- 
Kosterlitz-Thouless phase transition |Ber72l IKT73| : 
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Figure 18: The phase diagram of the 0{n) model on the hexagonal lattice. 



Conjecture 8.9. Let Xc{n) = l/\/2 + \/2 - n. For x < Xc{n) (resp. x > Xc{n)) the 
probability that two points are on the same loop decays exponentially fast (as a power 
law). 

The conjecture was rigorously established for two cases only. When n= 1, the critical 
value is related to the critical temperature of the Ising model. When n = 0, it was recently 
proved in jPCSlOj that ^2 + \/2 is the connective constant of the hexagonal lattice. 

It turns out that the model exhibits one critical behavior at Xc{n) and another on 
the interval (xc{n), +00)^ corresponding to dilute and dense phases (when in the limit 
the loops are simple and non-simple respectively), see Fig. 18 In addition to this, the 
two critical regimes are expected to be conformally invariant. 

Exactly as in the case of FK percolation, the definition of the spin fermionic observ- 
able can be extended. For a discrete domain with two points on the boundary a and 
6, the parafermionic observable is defined on middle of edges by 



F{z) 



-(TiW~^.{a,z) edges in loops in uj 

Q^-(TiWy(a,b)j^^ edges in ^fii^ loops in lu 



(8.4) 



where £{a,z) is the set of configurations of loops with one interface from a to z. One 
can easily prove that the observable satisfies local relations at the (conjectured) critical 
value if cr is chosen carefully. 

Proposition 8.10. If x = Xc - n, let F be the parafermionic observable 

with spin a = (7{n) = 1 - ^ arccos(-n/2); then 



{p-v)F{p) + {q-v)F{q) + {r-v)F{r) = 
where p, q and r are the three mid- edges adjacent to a vertex v. 



(8.5) 



This relation can be seen as a discrete version of the Cauchy-Riemann equation on 
the triangular lattice. Once again, the relations do not determine the observable for 
general n. Nonetheless, if the family of observables is precompact, then the limit should 
be holomorphic and it is natural to conjecture the following: 
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Conjecture 8.11. Let n e [0,2] and {Cl,a,b) be a simply connected domain with two 
points on the boundary. For x = Xc{n), 



where (7=1- ^ arccos(-n/2), Fs is the observable in the discrete domain with spin a 
and Tp is any conformal map from f2 to the upper half-plane sending a to oo and b to 0. 

A conjecture on the scaling limit for the interface from a to & in the 0{n) model can 
also be deduced from these considerations: 

Conjecture 8.12. For n e [0,2) and xdn) = l/\/2 + \/2 - n, as the mesh size goes to 
zero, the law of 0{n) interfaces converges to the chordal Schramm- Loewner Evolution 
with parameter n = 47r/(27r - arccos(-n/2)). 



This conjecture is only proved in the case n = 1 (Theorem 2.10). The other cases 
are open. The case n = is especially interesting since it corresponds to self-avoiding 
walks. Proving the conjecture in this case would pave the way to the computation of 
many quantities, including the mean-square displacement exponent; see |LSW04b] for 
further details on this problem. 

The phase x < Xc{n) is subcritical and not conformally invariant (the interface con- 
verges to the shortest curve between a and b for the Euclidean distance). The crit- 
ical phase X e (a;c(n),oo) should be conformally invariant, and universality is pre- 
dicted: the interfaces are expected to converge to the same SLE. The edge-weight 
Xc{n) = l/\/2 - \/2 - n, which appears in Nienhuis's works |Nie82l INie84| . seems to 
play a specific role in this phase. Interestingly, it is possible to define a parafermionic 
observable at Xc{n) with a spin (7{n) other than (7{n): 

Proposition 8.13. If x = Xc{n), let F be the parafermionic observable with spin a = 
(j{n) = -\ - ^ arccos(-n/2); then 

{p-v)F{p) + {q-v)F{q) + {r-v)F{r) = (8.7) 

where p, q and r are the three mid-edges adjacent to a vertex v. 

A convergence statement corresponding to Conjecture |8.11| for the observable with 
spin a enables to predict the value of k for Xc{n), and thus for every x > Xc{n) thanks 
to universality. 

Conjecture 8.14. For n 6 [0,2) and x e (xc(n),oo), as the lattice step goes to zero, 
the law of 0{ri) interfaces converges to the chordal Schramm- Loewner Evolution with 
parameter k = 47r/ arccos(-n/2) . 

The case n = 1 corresponds to the subcritical high-temperature expansion of the Ising 
model on the hexagonal lattice, which also corresponds to the supercritical Ising model 
on the triangular lattice via Kramers- Wannier duality. The interfaces should converge 
to SLE(6). In the case n = 0, the scaling limit should be SLE(8), which is space-filling. 
For both (slightly different) model is known to converge to the corresponding 

SLE (site percolation on the triangular lattice for SLE(6), and the perimeter curve of 
the uniform spanning tree for SLE(8)). Yet, the known proofs do not extend to this 
context. Proving that the whole critical phase {xc{n), oo) has the same scaling limit 
would be an important example of universality (not on the graph, but on the parameter 
this time). 

The two previous sections presented a program to prove convergence of discrete 
curves towards the Schramm-Loewner Evolution. It was based on discrete martingales 
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Figure 19: Different possible plaquettes with their associated weights. 



converging to continuous SLE martingales. One can study directly SLE martingales {i.e. 
with respect to cr(7[0, t])). In particular, g[{z)°''lgt{z) - Wt]^ is a martingale for SLE(k) 
where k, = 4(a - P)/[f3{f3 - 1)]. All the limits in these notes are of the previous forms, 
see e.g. Proposition |6.7[ Therefore, the parafermionic observables are discretizations of 
very simple SLE martingales. 

Question 8.15. Can new preholomorphic observables be found by looking at discretiza- 
tions of more complicated SLE martingales ? 

Conversely, in [SS05| . the harmonic explorer is constructed in such a way that a 
natural discretization of a SLE(4) martingale is a martingale of the discrete curve. This 
fact implied the convergence of the harmonic explorer to SLE(4). 

Question 8.16. Can this reverse engineering be done for other values of k in order to 
find discrete models converging to SLE? 



8.5 Discrete observables in other models 

The study can be generalized to a variety of lattice models, see the work of Cardy, Ikhlef, 
Riva, Rajabpour [I C091 IRC07I IRC06| . Unfortunately, the observable is only partially 
preholomorphic (satisfying only some of the Cauchy-Riemann equations) except for the 
Ising case. Interestingly, weights for which there exists a half-holomorphic observable 
which is not degenerate in the scaling limit always correspond to weights for which the 
famous Yang-Baxter equality holds. 

Question 8.17. The approach to two-dimensional integrable models described here is 
in several aspects similar to the older approaches based on the Yang-Baxter relations 
jBaxS^. Can one find a direct link between the two approaches? 

Let us give the example of the 0{n) model on the square lattice. We refer to ITOOS] 
for a complete study of the following. 

It is tempting to extend the definition of 0{n) models to the square lattice in order to 
obtain a family of models containing self-avoiding walks on Z'^ and the high-temperature 
expansion of the Ising model. Nevertheless, difficulties arise when dealing with 0(n) 
models on non-trivalent graphs. Indeed, the indeterminacy when counting intersecting 
loops prevents us from defining the model as in the previous subsection. 

One can still define a model of loops on G c L by distinguishing between local 



configurations: faces of G* c L* are filled with one of the nine plaquettes in Fig. 19 A 
weight py is associated to every face v e G* depending on the type of the face (meaning its 
plaquette). The probability of a configuration is then proportional to n'^ Ilveh* Pv- 
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Remark 8.18. The case ui = U2 = v = x, t = 1 and wi = W2 = n = corresponds to 
vertex self-avoiding walks on the square lattice. The case ui = U2 = v = \/w^ = \/w2 = x 
and n = t = 1 corresponds to the high-temperature expansion of the Ising model. The 
case t = Ui = U2 = V = 0, Wi = W2 = 1 and n > corresponds to the FK percolation at 
criticality with q = n. 



A parafermionic observable can also be defined on the medial lattice: 



F{z) 



y N (.-i<^W^(a,z) # loops i-r 



where £ corresponds to all the configurations of loops on the graph, and £{a,z) corre- 
sponds to configurations with loops and one interface from a to z. 

One can then look for a local relation for F around a vertex v, which would be a 
discrete analogue of the Cauchy-Riemann equation: 

F{N)-F{S) = i[F{E) - F{W)l (8.9) 

An additional geometric degree of freedom can be added: the lattice can be stretched, 
meaning that each rhombus is not a square anymore, but a rhombus with inside angle 
a. 

As in the case of FK percolations and spin Ising, one can associate configurations by 



pairs, and try to check (8.9) for each of these pairs, thus leading to a certain number of 
complex equations. We possess degrees of freedom in the choice of the weights of the 
model, of the spin a and of the geometric parameter a. Very generally, one can thus try 
to solve the linear system and look for solutions. This leads to the following discussion: 



Case V = and n = 1: There exists a non-trivial solution for every spin a, which 
is in bijection with a so-called six-vertex model in the disordered phase. The height 
function associated with this model should converge to the Gaussian free field. This 
is an example of a model for which interfaces cannot converge to SLE (in f IC09| : it is 
conjectured that the limit is described by SLE(4,p)). 



Case V = and n i= 1: There exist unique weights associated to an observable with spin 
-1. This solution is in bijection with the FK percolation at criticality with ^^q = n + 1. 
Nevertheless, physical arguments tend to show that the observable with this spin should 
have a trivial scaling limit. It would not provide any information on the scaling limit of 
the model itself; see |IC09| for additional details. 



Case V i= 0: Fix n. There exists a solution for o" = |^ - | where rj e [-7r,7r] satisfies 
-■^ = cos2?7. Note that there are a priori four possible choices for a. In general the 
following weights can be found: 



t 


= -sin(2(/>- 


3r]/2) + sin(5r;/2) 


Ui -- 


= -2 sin(7y) 


cos(377/2 - (j)) 


U2 -- 


= -2 sin(77) 


sin((/)) 


V 


= -2 sin((/)) 


cos(3r7/2 - (j)) 


Wi -- 


= -2 sm{(j)- 


- ij) cos(377/2 - (j)) 


W2 = 


-- 2 cos(?7/2 


- (j)) sin((/)) 



sin(3r//2) + sin(ry/2) 



where = {1 + a)a. We now interpret these results: 
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When 7] e [O.tt], the scaling limit lias been argued to be described by a Coulomb gas 
with a coupling constant 2r]/Tr. In other words, the scaling limit should be the same as 
the corresponding 0(n) model on the hexagonal lattice. In particular, interfaces should 
converge to the corresponding Schramm-Loewner Evolution. 

When r] € [-7r,0], the scaling limit curve cannot be described by SLE, and it pro- 
vides yet another example of a two-dimensional model for which the scaling limit is not 
described via SLE. 
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